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AN  APPROXIMATION  THEORY  FOR 
NONLINEAR  PARTIAL  DIFFERENTIAL  EQUATIONS 
WITH  APPLICATIONS  TO  IDENTIFICATION  AND  CONTROL* 

H.  T.  BANKSt  and  K.  KUNISCHt 


Abstract  Approximation  results  from  linear  searifrotip  theory  are  uwd  to  develop  a  genera!  framework 
for  cooverjenoe  of  approximation  tchemes  in  parameter  estimation  and  optima]  control  problems  for1 
nonlinear  partial  differential  equations.  These  ideas  are  used  to  establish  theoretical  convergence  results 
for  parameter  identification  using  modal  (eigenfunction)  approximation  techniques.  Results  from  numerical 
investigations  of  these  schemes  for  both  hyperbolic  and  parabolic  systems  are  given. 
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1.  Introduction.  When  modeling  real-world  phenomena  one  often  encounters  a 
situation  where  a  priori  knowledge  leads  one  to  conjecture  a  certain  type  of  model 
equation  containing  parameters  which  are  unknown.  In  this  paper  we  are  primarily 
concerned  with  technique*  for  recovery  of  these  unknown  quantities  from  given  data. 
In  99  2  and  3  we  present  a  quite  general  framework  for  approximation  schemes  for 
abstract  nonlinear  Cauchy  problems.  These  approximation  results  are  subsequently 
applied  to  modal  techniques  for  identification  and  control  problems  in  95  4  and  5, 
respectively.  A  summary  of  some  of  our  numerical  experience  with  parameter  estima¬ 
tion  problems  using  these  techniques  is  given  in  5  6.  The  examples  here  were  chosen 
so  as  to  illustrate  the  feasibility  and  effectiveness  of  the  method  and  to  investigate 
possible  inherent  difficulties.  We  are  quite  confident  that  the  ideas  outlined  here  will 
he  applicable  in  a  variety  of  research  areas  where  mathematical  models  for  the 
phenomena  under  study  are  used.  In  a  forthcoming  monograph  we  shall  discuss  in 
more  detail  identification  problems  that  arise  in  several  areas  of  applications  [35] 
including  seismology  [31  [101  [181  reservoir  engineering  [111  [17],  [381  glaciology 
[161  physics  [371  biology  [41  [51  [291  £34]  and  large  space  structures  technology. 
While  our  treatment  here  is  restricted  to  constant  unknown  parameters,  the  theoretical 
ideas  extend  in  large  part  to  problems  with  unknown  function  parameters.  Indeed, 
we  are  currently  applying  some  of  our  techniques  to  specific  problems  from  the  areas 
mentioned  above;  in  some  cases  these  efforts  involve  identification  of  functions. 

In  this  paper  die  general  approximation  results  are  used  to  carefully  discuss  modal 
approximation  schemes  for  certain  classes  of  partial  differential  equations  (see  59  4 
and  5).  Such  schemes  for  specific  identification  and  control  problems  are,  of  course, 
not  new.  Many  discussions  in  the  literature,  however,  are  in  the  context  of  very  specific 
examples  and  frequently  no  convergence  proofs  or  evidence  of  numerical  studies  are 
supplied.  Modal  approximations  have  many  advantages,  including:  they  are  readily 
discussed  and  understood  in  terms  of  classical  spectral  results;  they  are  familiar  to 
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and  readily  implemented  by  practicing  engineers,  and  they  give  rise  to  a  simple 
algebraic  structure  for  the  approximating  ordinary  differential  equations.  However, 
modal  approximations  do  have  some  shortcomings,  of  which  we  mention  several. 
First,  in  many  practical  problems  it  is  very  difficult  to  calculate  the  true  natural  modes. 
Secondly,  for  certain  parabolic  partial  differential  equations  modal  approximations 
by  their  very  nature  lead  to  stiff  systems  of  approximating  ordinary  differential 
equations.  Finally,  one  can  encounter  lack  of  “numerical  identiffability’’  (i.e.,  the 
identification  problems  for  the  approximating  ordinary  differential  equations  yield 
parameter  estimates  that  converge  to  different  values  for  different  sets  of  initial 
estimates)  regardless  of  the  well-posedness  of  the  parameter  estimation  problem  for 
the  original  partial  differential  equation  model.  With  respect  to  the  first  difficulty 
pointed  out  here,  we  refer  to  Example  4.4,  below,  where  we  explain  a  “modal” 
approximation  scheme  for  an  identification  problem  which  does  not  employ  the  natural 
modes  of  the  system.  For  one  solution  of  the  latter  problems,  our  experience  indicates 
that  for  certain  classes  of  parabolic  problems  spline -based  approximation  schemes 
can  be  more  efficient.  Details  on  this  aspect  of  spline  methods,  along  with  a  number 
of  other  features  of  these  techniques,  will  be  given  in  a  separate  manuscript  currently 
in  preparation. 

The  parameter  identification  and  estimation  problem  has  received  a  great  amount 
of  attention  in  the  engineering  literature  and  we  refer  to  [1],  [23],  [31],  [32]  for  review 
articles.  In  the  future  monograph  alluded  to  above,  we  shall  survey  the  research  efforts 
from  the  engineering  as  well  as  from  tile  mathematical  literature.  Much  of  the 
mathematical  literature  is  concerned  with  the  problem  of  identifiability,  which,  loosely 
speaking,  is  defined  as  the  problem  of  injectivity  of  the  map  from  the  set  of  parameters 
to  the  set  of  outputs.  Although  this  is  a  very  important  theoretical  and  practical 
question,  it  will  not  be  a  part  of  the  discussion  of  the  present  paper. 

We  point  out  one  important  technical  aspect  that  will  become  clearer  in  Examples 
4.1  and  4.4  below.  In  general,  the  eigenfunctions  of  the  model  equation  will  depend 
on  the  parameters  that  are  to  be  identified.  For  modal  approximation  schemes  this  is 
an  extremely  undesirable  feature  from  the  point  of  view  of  implementation,  since  in 
practical  examples  the  representation  of  the  operators  in  the  approximating  equations 
will  involve  a  matrix  of  inner  products  of  the  eigenfunctions.  It  is,  of  course,  desirable 
to  have  this  matrix  independent  of  the  unknown  parameters  to  avoid  excessive 
numerical  integrations  when  performing  iterative  seances  on  these  parameters. 

Our  focus  in  this  paper  is  on  the  development  of  semkiiscrete  approximation 
schemes  for  parameter  identification  and  control  problems  which  result  in  approximat¬ 
ing  problems  governed  by  ordinary  differential  equations.  Of  course,  full  discretgation 
methods  (dkcretixstion  in  time  as  well  as  spatial  coordinates,  resulting  in  problems 
governed  by  difference  equations)  are  of  great  importance  and  our  investigations  of 
a  related  theoretical  framework,  as  well  as  detailed  schemes  for  such  an  approach, 
will  be  reported  elsewhere. 

In  summary,  the  emphasis  of  our  presentation  is  twofold.  First,  we  present  a 
genera  toeoreocii  crane  worn,  wits  munown  ptruDetef-oepenoeBt  tpeoet,  wmcn 
can  be  used  to  tram  many  types  of  problems  (including  estimation  of  function  space 
parameters)  and  approximation  schemas  (sat  the  remarki  in  1 7  below).  As  a  concrete 
example  of  the  nee  of  tide  framework,  we  give  a  detailed  treatment  of  “modal” 
approximation  schemas,  thereby  putting  on  a  found  theoretical  foundation  a  dam  of 
methods  that  have  beau  mad  hi  an  ad  hoc  wty  by  tcienttofs  and  engineers  for  some  time. 

The  notation  need  throughout  the  paper  to  write  standard.  We  employ  the  usual 
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denote  norms  of  elements,  as  well  as  those  of  operators.  Only  in  cases  where  confusion 
may  arise  will  we  use  subscripts  to  distinguish  norms  in  various  spaces. 

2,  The  abstract  tdanttlraHoa  problem  and  Us  apprnriastlnn  We  consider  the 
abstract  semilinear  Cauchy  problem 

u(t)=A(q)u(t)+F(q,  t,  u(t)).  t>0, 

(2.1) 

u(0)  =  uo(q), 

where  for  each  qeQc:Rk,A(q)ulb6  infinitesimal  generator  of  a  linear  Co-  semigroup 
{?(f;q)}IB0  on  a  real  Hilbert  space  X(q)  with  inner  product  (•,•},  and  norm  |-|, 
(denoted  sometnnes  below  by  X,  <•,  •)  and  |>|,  respectively,  when  no  loss  of  clarity 
results  from  suppression  of  the  q).  We  shall,  throughout  our  discussions,  employ  the 
concept  of  mild  solutions,  so  that  t-*u(t;q)  is  called  a  solution  of  (2.1)  if  it  satisfies 

(2.2)  u(f,q)~T(.t;q)ual(q)+^  T(f-s\q)F(q,s,u(s;q))ds. 

We  note  that  for  solutions  u  we  have  t-*u(t;q)  continuous.  The  conditions  that  we 
impose  on  F  below  will  guarantee  existence  and  uniqueness  of  mild  solutions  u  of 
(2.1)  on  any  given  finite  interval  [0, 1”].  We  shall  in  certain  specific  instances  below, 
be  required  to  discum  briefly  die  relationship  between  mild  and  strong  (in  a  classical 
almost-everywhere  sense)  solutions  of  (2.1),  but  for  more  general  results  are  refer  the 
reader  to  [28]. 

Throughout  our  presentation  we  shall  assume  that  X(q)  is  a  function  space  of 
/?" -valued  “functions’*  (possibly  one  of  the  usual  Lebesgue  spaces  of  equivalence 
dames  of  functions)  defined  on  the  fixed  interval  [0, 1];  consequently,  we  shall  also 
use  the  notation  u(t,x;  q)  or  u(t,  ■ ;  q)  when  discussing  solutions  at  (2.1). 

While  we  shall  also  discum  control-theoretic  applications,  much  of  our  attention 
arill  be  directed  towards  the  problem  of  identifying  the  parameter  q  in  (2.1)  from 
observations  of  the  system.  Specifically,  we  assume  tint  (2.1)  models  some  physical, 
biological,  economic,  etc.,  system  for  which  output  measurements  j)  are  available. 
These  measurements  may  be  available  in  the  form  of  continuous  data  f(t),  Q&t&T, 
or  discrete  data  fit,),  0Sf><-  •  •  <t,sT.  We  then  seek  to  find  a  “best”  value  for  q 
in  Q  by  minimizing  an  appropriately  defined  fit-to-data  criterion.  To  be  specific  in 
oar  formulation  fame,  we  shall  assume  discrete  time  observations  with  values  f(t,)  in 
an  observation  space  9.  All  of  the  results  of  this  paper  are  easily  extended  to  the  case 
of  identification  problems  where  one  hm  continuous  time  data,  but  we  shall  not  pursue 
such  problems  here.  Assuming,  then,  that  a  criterion  function /:<?x C(0,  T-,X(q)) x 
n^.i  9+R1  is  defined,  we  formally  state  the  identification  problem: 

(ID)  Given  observations  f  ~{?0s)K-i,  minimize  f(q,  u(-;q),  f)  over  qeQ  subject 
to  u( • ;  q)  satisfying  (2.2). 

Several  traditional  choices  of  fit-to-data  criteria  are  included  in  oar  formulation; 
namely,  we  may  consider  either  integral  or  pbiutwise  (in  a  spatial  sense)  evaluation 
least-squares  sums  in  the  above  formulation.  In  the  case  al  integral  evaluation  we  are 
given  measurements  ?(4)eLs(0,l‘,jr)  where  *£*  and  an  output  map 
Y(t,x,q):RH-*R''  on  the  “state”  u(t,x;q).  The  observation  space  is  given  by  •  « 
1^(0, 1;H")  and  the  criterion  is  defined  by 

(2.3)  /<« [,|f(h,s)-r^*,«)h(ls*tq)!a<fc. 

imt  Jo 
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We  assume  that  Y  is  continuous  is  q  and  sufficiently  regular  in  x  so  that  x-* 
Y((i>x,q)u(thx;q)  k  in L2( 0,  UR”).  For  the  choice  of  poiotwise  or  spatially  discrete 
measurements,  we  assume  that  we  have  observations  y0i)e&  »nj-i  R",  correspond¬ 
ing  to  measurements  of  the  output  at  points  in  [0,  l]at  time  //.These  observations 

represent  measurements  for  C(ihq)((thq)  where  ((th  q)  =  coi  (u(t,,  xt;  q),  •  •  ■ , 
u(th  x,;  q))  and  C(t,,  q)  is  an  (*■/)  x  («/)- matrix  depending  continuously  on  q  for  each 
fixed  t,.  The  associated  fit-to-data  criterion  is  then  defined  by 

(2.4)  /(«.«(•  ;«),*)“  Z  [m-C^qmuqtf. 

i 

The  output  maps  Y  and  C  introduced  hi  (2.3)  and  (2.4)  are  necessitated  by  the 
fact  that  often  in  practice  one  can  observe  only  some  components  (say  v)  of  the 
n- dimensional  vectors  u({,  x ;  q),  and  that  these  observations  may  depend  on  the  time 
at  which  they  are  made.  We  further  note  that  the  point  evaluations  at  xt  used  in 
defining  f(fh  q)  above  may  be  meaningless  without  additional  restrictions  on  the  state 
space,  the  initial  data  and/or  die  right  side  of  the  equation  in  (2.1).  A  more  detailed 
discuasion  of  the  problems  arising  from  use  of  criteria  such  as  (2.4)  when  dealing  with 
mild  solutions  will  be  given  in  the  context  of  Example  4.3  below. 

We  turn  next  to  formulating  a  sequence  (IDV)  of  approximating  problems  on 
Hilbert  spaces  Xs (q)  for  our  original  identification  problem  (ID).  These  problems 
involve  “states”  governed  by  ordinary  differential  equations  and  are  (in  the  specific 
instances  we  shall  propose)  tractable  using  standard  numerical  procedures.  We  state 
first  a  series  of  hypotheses  and  definitions  that  will  be  needed  at  various  points  in  the 
sequel. 

(HI)  For  each  JV-1,2,  •  •  •  ,XN(q)  is  a  closed  linear  subspace  of  X(q),  endowed 
with  the  X(q)  topology. 

(H2)  The  spaces  X(q),  qeQcRk,  are  set-theoretkally  equal  and  uniformly 
topologically  isomorphic  so  that  there  exists  a  constant  JTSl  such  that  Iv^S 
3T|e|,  for  aO  v  eX -X(q)  and  q,  4 € Q. 

(H3)  For  each  q  a  Q,  A(q)  generates  a  Unear  Co-semigroup  T(t;  q)  on  X(q). 

(H4)  The  set  Q  k  a  compact  subset  of  R \ 

(H5)  (i)  For  each  <?€<?,  let  P* (q):X{q)-*XN{q)  denote  the  canonical  orthogonal 
projections  along  and  let  A"(q):X(q)-*  XN(q)be  defined  by  A*(q)- 

Ps(q)A{q)PN (q).  For  each  N,  let  Aft(q)  generate  a  linear  Co-semigroup  on 
JT(q)  denoted  by  r*(f ;  q). 

(ii)  For  each  N,  there  exist  constants  A  mA(N)  and  <5  »<5(A0,  independent 
of  q,  such  that  |T*(f ;  q)j  sA  e *. 

(H6)  (i)  For  each  continuous  function  u:  (0,  T]-*X  ~X(q)  (see  (H2)),  the  map 
t  -*F(q,  t,  u(f ))  is  measurable. 

(ii)  For  each  constant  M>0,  there  exists  a  function  ktmkt(Af)  in  Lj( 0,  T) 
such  that  for  any  q,4*0  we  have 

lF(q,  U  «i) -F(q,  t,  Ma)|s  S  *j(»)|«i  -  uiU 

(Or  afi  Hi,  Ha«JT  with 

(Hi)  There  exists  a  function  kj  in  Li(0,  T)  such  that 

\Fiq,  k  e)lsffi*a  (fMMs  +  1} 

for  all  ve X, q,4* Q. 

(hr)  For  each  (l,  p)o[0,  T]xX,  the  atap  q  -*F{q,  t,v)  k  continuous.  (Again, 
saadar  (H2)  recall  that  aO  tha  sats  X  -  *Cq)  are  the  same.) 


APPROXIMATION  OF  NONLINEAR  POES 


819 


(H7)  The  projections  P* (q)\X(q)-*XN (q)  are  such  that  for  any  sequence  {qN}  in 
Q  satisfying  -*<f  e  Q,  one  has  ji’7''  (qN)z  -z),*  -►()  as  Af-*aofor  all  r  eX(Q). 
(H8>  For  each  convergent  sequence  qN  -*<?  in  Q,  there  are  constants  M,  m  such 
that  lrN(r;q")U &M  e~,  |T(r; flj,  S  Mem  uniformly  in  N  - 1, 2,  •  •  • . 

(H9)  For  each  convergent  sequence  q*  -*  q  in  Q,  one  has  for  z  eX(4),  |rN(r ;  qN)z  - 
T(f,  4)*^  -*0  as  N  •*  uniformly  in  t  e  [0,  T\ 

The  assumption  (H2)  mil  be  taken  as  a  standing  hypothesis  for  the  remainder  of 
our  discussions.  For  the  approximating  schemes  we  develop  below,  consistency  will 
follow  from  (H7)  while  (H8)  is  a  statement  of  stability.  As  we  shall  see,  convergence 
of  the  schemes  (which  is  (H9))  will  follow  from  (H7),  (H8)  and  the  Trotter-Kato 
theorem. 

Remark  2.1.  It  suffices,  under  the  standing  assumption  (H2),  that  the  following 
condition  hold  in  place  of  (H6)(ii):  For  some  fixed  q*eQ  we  have  that  for  each 
Af >0  there  is  a  function  kx  such  that  tor  all  qeO  the  relation  |F(q,r,ui)- 
F(qtt ,  uz)(,*  S  k  i(()|ut  -  Mat,*  feu  all  uu  u2eX  with  |  u/^SAf.  Indeed,  it  is  easily  seen 
that  this  condition,  along  with  (H2),  implies  (H6)(ii).  Similarly,  we  can  in  the  presence 
of  (H2)  equivalently  postulate  in  place  of  (H6)(iii)  the  conditions:  For  some  fixed 
q*eO  there  exists  a  function  k2  such  that  tF(q,  t,  S  *j(r  MM*- + 1}  for  all  veX, 

qeQ.  We  further  note  that  existence  of  a  function  kiZL2(Q,  T)  such  that  |F(q,  r,  0)|„-  £ 
kilt)  for  q,qeQ,  a  statement  of  the  inequality  of  (H6)(iii)  holding  only  for  |u|« 
sufficiently  large  (i.e.t  affine  growth  at  ao),  along  with  (H6)(ii),  are  sufficient  to  imply 
(H6*Bi). 

While  the  complete  role  played  by  the  various  hypotheses  in  our  development 
will  be  dearer  after  our  presentation,  a  few  explanatory  comments  here  might  be 
helpful  te  readers.  First,  the  desirability  of  the  generality  of  allowing  the  underlying 
Hilbert  space  X  for  (2.1)  to  depend  on  q  in  such  a  way  that  (H2)  obtains  will  not  be 
apparent  from  the  examples  discussed  here.  (Rather,  one  must  for  this  consider  certain 
parabolic  problems— see  the  comments  in  9  7.)  However,  in  light  of  (H2>  as  a  standing 
assumption,  we  ate  justified  in  suppressing  the  canonical  isomorphism  JN  :X($)-» 
X(qN)  in  writing  \Pff< («*')*  -r[^v-*0  in  (H7)  rather  than  the  technically  correct  state¬ 
ment  \PN(qN)J^,z Similar  observations  are  pertinent  for  the  statement 
of  (H9)  as  well  as  in  numerous  other  places  in  our  presentation  where  we  suppress 
the  notation. 

Condition  (H4),  while  seemingly  stringent,  is  an  assumption  often  valid  in  practical 
problems  where  our  theory  might  be  useful.  Since  under  (H3)  A(q)  is  closed,  it  follows 
from  the  closed  graph  theorem  that  AN(q)  of  (H5)(i)  is,  in  fact,  bounded  and  hence 
(H5Xi)  follows  immediately  from  (H3).  It  should  be  recognized  that  the  form  of  the 
approximating  operators  defined  in  (H5)  is  a  classical  one  (e.g.,  see  [33,  p.  3695  which 
has  aho  recently  been  employed  in  the  development  of  spline  approximation  tech¬ 
niques  for  delay  differential  equations  [6].  The  definition  at  AN(q)  involves  the  implicit 
assumption  that  Xs (q)  <=  Dom  (A  (<?));  since  our  goal  here  is  the  rigorous  formulation 
of  modal  approximation  schemes  for  (ID),  this  restriction  poses  no  difficulties. 
However,  it  does  prevent  a  straightforward  inclusion  of  low-order  finite-element 
metitods  for  higher-order  partial  differential  equations  in  our  approximation 
framemort. 

Hypothesis  (H7)  is  a  common  requirement  (e.g.,  see  [24],  [305  in  approximation 
dsnory,demanding  that  the  sequence  X’N  of  subspaces  actually  approximate  the  original 
state  space  X.  Ftosfty,  916)  is  comprised  of  conditions  on  the  nonHnearities  in  (2.1) 
Uiii  ire  funGwony  hiiiiii  ro  incrooo  muiy  uuwnuii  pruoicixiB  oc  prscocu  lmpor- 

tance  but  yet  are  strong  enough  to  guarantee  global  existence  of  solutions  of  (2.1)  on 
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fixed  finite  intervals.  As  the  knowledgeable  reader  might  expect,  these  conditions  can 
be  replaced  by  alternate  and/or  weaker,  hypotheses,  but  only,  in  general,  at  the  cost 
of  xtriirinrml  tedium  in  the  proofs  below.  We  have  tried  to  compromise  between 
strong  conditions  that  are  easily  stated  and  employed  in  the  proofs  and  ones  that  are 
as  general  (and  weak)  as  possible.  Further  comments  on  this  matter  will  be  made  in  i  3. 

Before  defining  the  approximating  equations  for  (2.1),  we  define  the  projection 
of  the  nonlinearity  F  onto  XN  by  FN(q,t,o)»PN(q)F(q,t,v)  for  each  (q,t,  c)e 
Q  x [0,  T\y-X.  The  approximating  family  of  equations  is  then  given  by 

v(t)=AN(q)v(t)+FN(q,t,v(t)),  t>  0, 

(2‘5)  v(0)=PN(q)uo(q). 

Assuming  existence  of  (mild)  solutions  to  (2.5)  (tins  will  be  established  below),  we 
denote  (for  a  given  q)  these  solutions  by  u*(f)  or  alternatively  uM(t;q)  or  uN(t,  x ;  q), 
depending  on  the  context.  We  then  define  the  approximate  identification  problems 
(ID*)  by: 

(ID*)  Given  observations  y  »{y(fi)K-i  and  a  fit-to-data  criterion  /,  minimize 
JN(q)»/(q,  uN(-:q),  y )  over  qeQ  subject  to  u*(-;<?)  satisfying  (2.5). 

We  note  thtf  if  (in  addition  to  (HI ))XN(q)  is  finite-dimensional,  then  (2.5)  can  be 
equivalently  interpreted  in  the  strong  sense  and  (ID*)  then  becomes  an  optimization 
problem  constrained  by  finite-dimensional  ordinary  differentia]  equations. 

In  our  discussions  below,  we  shall  denote  by  <f*  any  solution  of  (ED*)  so  that  it 
follows  by  definition  that  /*(4*)5/*(q)  for  all  q  e  Q. 

Proposition  2.1.  Assume  that  (H2),  (H3)  and  (H6)  obtain.  Then  for  each  qeO 
there  exists  a  unique  (mild)  solution  «(•  ;q)eC(0,  T;X(q))  of  (2.1).  If,  in  addition, 
(H5)  holds,  then  exists,  for  each  N ■  1,2,  •  •  • ,  a  unique  (mild)  solution  u  (• ; q)e 
C(0,T;X(q))  of  (2.5). 

Proof.  The  proofs  are  completely  standard  and  we  only  sketch  the  ideas  tor  (2.1). 
Uniqueness  follows  immediately  from  (H6)  and  an  application  of  Gronwall’s 
inequality.  Existence  is  established  through  the  usual  Picard  iterate  techniques.  Define 
o°(f)  -  T(r;  q)u0(q)  and  for  /  - 1, 2,  •  •  • , 

(2.6)  o'(t)-r(f;q)siofa)+J[  T(t-s;q)F(q,s,v,~1(s))ds 

for  t €[0,  T).  Prom  (H3)  and  (H6)  it  is  easily  seen  that  the  iterates  v1  are  all  well 
defined  and  o':  [0,  T]-*X(q)  is  continuous.  Moreover,  {v‘)T-o  is  a  bounded  subset  of 
C(0,  T;X).  Employing  (H6Kii)  and  simple  inductive  arguments,  one  can  establish 
diet  {o'}  is  Cauchy  fat  C(0,  T;  X).  Panting  to  the  limit  in  (2.6),  one  obtains  the  desired 
results.  Existence  of  unique  solutions  of  (2.5)  is  argued  in  an  analogous  manner  by 
appealing  to  (H5)  for  appropriate  boundedness. 

Thsomm  2.1.  Assume  hypotheses  (H1MH6)  hold  and  let  /(*,•,?): 

C(0, T-,X)-*Rl  be  continuous.  Moreover,  suppose  q-*«ofa),  q-*F*(q)z  and  q-* 
TN(t,q)s,xvX,  are  continuous,  with  the  latter  uniformly  In  ta[0,Tl  Then:  (i)  There 
exists  for  each  # a  soiuOon  4"  of  (ID*)  and  the  sequence  M  )  possesses  a  convergent 

subsequence  4*k-*4-  (H)  If  we  further  assume  that,  for  any  sequence  fa'}  in  Q  with 
,  w*  have  |e'(r;q') -»0  «#/■*«  uniformly  bt  ta[0,  TJ,  then  4 1*  * 
solution  of  (JD). 
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Proof.  We  show  for  fixed  N  that  q  ■*JN(q)*J(q,  mn(-  ;  q),  y)  is  continuous  on  Q 
which,  in  the  light  of  (H4),  yields  (i).  First  note  that  uN  satisfies 

(2.7)  uN{r,q)  =  TN(t;q)PN(qMq)+  f  TN(t-s,q)PN(q)F(q,s,uN(s;q))ds 

Jo 

on  [0,  n  1°  view  of  (HS)  and  (H2)  we  find  that  there  exist  constants  M  and  Mi  such 
that  jr"^;  q%  Silf,  SMi  uniformly  tor  t  £  [0,  T\  q,  q  e  Q.  It  follows  from 

(H6)(iii)  and  (2.7)  that 

|«w(t;  q%  SAflPN(q)«0(q)|,  +  M  £  \PN{q)F{q,  s,  uN(s ;  q%,  ds 
SAf |«0(4)I+m|o  k2(s){\uN(s\ q)|,  + 1}  ds. 

Since  q-»Uo(q)  is  continuous,  this  implies  (via  Gronwall’s  inequality)  that  |uv(f;q)|, 
is  uniformly  bounded  for  (t,  q)  €  [0, 71]  x  Q.  This  in  turn  implies  (by  (H6)(ih))  that  the 
mapping  t -*  TN(t-j;q)PN(q)F(4,  s,  uN(s\4))  horn  [0,  T]  to  X  is  dominated  by  an  av 
integrabk  function  uniformly  in  q,q,tfeQ.  (This  will  permit  us  to  invoke,  below,  /.ijJSv 
the  usual  dominated  convergence  theorem.)  Assuming  that  q'-*q,  q‘,  tfeQ  are 
arbitrary,  we  obtain  the  following  estimates:  / 


|«w(r;4)-«K(r;q/)| 

;  4)Pn  (4)uc(4)~TN(t;  | 

+|Tw(r;  q')PM(l4M4)  -  T"(f,  q')P'V)«o(q)| 
+|TN(r;  <V,dl)«s(D  ~  TN(t;  JiP*  (q‘)uo(q‘)\ 


+|rw(r;  -  rK(r;  q,)P*V)«o(<l*)| 

+1  TN(r;  <?')P'V)«o(4)  -  TN(t;  q'lP'V  W)  | 

+  f  |{TN(f - j; 4) - TN(t-s; q‘)}PNW(q, s,  uN(s;q))\ds  /*•* 

*  I . JO  / 

+£  \TN(t-s;q‘)(PN(4)-PN(qlm4,  *,  uw(s;4))|dr  /  ?•* '** 


+jf  \T"(f  -  a;  q'^iq’mi  r,  h  w(r ;  4))-F(q‘,  s,  us(s ;  q'))}\  ds 


'o  r* 
2  / 


-  Pi(/) +p2(/)+Ps(/)+P*(/)  +Ps  (/) +P«(/), 

where  the  p(’s  are  defined  as  indicated  (p(  the  rth  term),  i  - 1,  •  •  • ,  6,  and  all 
are  |  •  l*.  We  then  have  by  hypothesis 

*(/) - Hrw«; 4)-r*(/; q')}i^(<n«o(4)h 0  as/-*®, 

uniformly  in  r  e[0, 1*].  Also  pa  and  ps-*0  by  tile  continuity  assumptions  an  Ps  and 
Ho  and  die  boundedness  of  TN(t;q').  Dominated  convergence  implies  that  p4-*0 
and  p5  ■*  0  as  /  -*  w.  Finally, 

Po(/)*A«fi|£  I F(4,  r,  u"(s-,4))-F(q‘,  s,  u"{,s,4))\ds 

+£  W, »,  H*(r ;  i))-F(^,  s,  uw(s;q,))|  dr} 
SMWjprC/H^  *»(r  )|h  w  (r ;  $)  -  «  "  (r ;  q')l  dr, 
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where  (h(j)-*0  as/-»oo  and  ki  depends  on  die  uniform  bounds  for  uN(t-,q),  uN(f>q  ) 
(see  (H6)(ii)).  Thus  we  fold 

|KK(f;4)-«l'r(*;g,)|£e/(f)+£  *i(r)|i«N(j;4)-K,'(j;$,)|d*, 


where  «;(<)-*  0,  uniformly  in  t  e  [0,  T],  as  /-*  oo.  Applying  Gron wall’s  inequality,  we 
have  uniform  (in  t)  continuity  otq-*uN(f,q)oaQ  which  implies  die  desired  continuity 
of/". 

Turning  to  (u)  and  letting  (4*'}  be  a  convergent  subsequence  with  $N‘  -*q,  we 
first  observe  that/JV'(4N0  SJNl(q)  for  all  q  e  Q.  By  hypothesis,  uN'(t;  QN<)-*u(f,  q)  and 
furthermore  «*'(*;  q)-*u(t;  q)  for  each  qsQ,  with  convergence  in  both  cases  uniform 
in  t  on  [0,  T\  This  indies  fN,(4N,)-*J(4, u(- ,4),  f)  and  fN'(q)-*J(q, u(  ;q), y)  as 
/-*oo  and  hence  from  the  above  inequality  we  obtain,  by  passing  to  the  limit, 
/(&«(•;  4),  f )  S f(q,  «(■;  q),  f)4or  utt  qeQ.  Thus  4  is  a  solution  of  (ID)  and  Theorem 
2.1  is  estabhrited. 


Remark  2.2.  If  e  -*J(q,  v,  #)  as  a  mapping  on  C( 0,  T ;  X)  actually  depends  only 
on  a  finite  number  of  values  v(tt),  t,  e  [0,  T\  as,  for  example,  in  (2.3)  or  (2.4),  then 
the  hypotheses  of  Theorem  2.1  involving  uniformity  in  f  can  be  relaxed  to  statements 
^holding  only  for  each  fixed  t  e  [0,  T].  The  above  arguments  remain  unchanged  except 
at  the  uniform  (in  /)  convergence  remarks  are  replaced  with  pointwise  convergence 
ements  (see  especially  the  term  pi(f)  in  dm  proof). 

^We  conclude  this  section  with  a  brief  explanation  of  how  (2.S)  (or  (2.7))  is  to  be 
.  actual  computations.  We  adopt  notation  very  similar  to  that  found  in  [6]  in 
elopment  of  spline  methods  for  delay  systems.  We  assume  that  XN  is  finite- 
and  choose  a  basis  independent  of  q  (recall  that  (H2)  is  a  standing 

p)by 

0s  m  (fiu  •  ’  •  •  P<nm) 

pdfaa XN(q).  From  (2.7)  under  (H5)  we  see  that  the  solution  uN  of  (2.5) 
for  aB  t  and  hence  there  edits  a  representation  uN(t\q)^ 
w  J^TwiA  w"(t ;  q )  -  col  (w?(t;q),  •  •  • ,  wJW)(f;  q))eR*°°.  We  let  [AN(q)] 
and  E r(<9.  t,  w  )]  denote  the  matrix  and  coordinate  representations  relative  to  0N  of 
AN(q)  and  P*{q)F{q,  t,  #*),  respectively.  The  coordinate  representation  of  (2.5)  is 
then  given  by 


^  <l)~[XN{q)]wN(ti  q)+[FN(q,  t,  w"(f,  q))\  r>0, 

w*(0;q)«r  , 

*  vdtere  ys  is  dstned  through  -/f^w"©;*).  For  any  zeXN  the 

associated  coordinate  vector  a"**4™  in  P"(q)e -0"a"  k  determined  uniquely 
by  the  comBtioa  (F"(q)X“t)lX",  or  equivalently 
Thus  we  have 

(2.9)  as-(Q"riRNz 

where  Q"  k  tht  4^nxd(ff)  matrix  with  elements  <0?,0?),  and  (R*r),  - <0?,  r), 
for  /  - 1, 2,  •  •  • ,  d(A0.  Arguing  euately  as  In  [6,  pp.  508-509],  we  therefore  find 

<2.io)  J  €**&»*«*>“***, 

where  X"  k  the  4W*dW)  matrts  with  etemeMs  K?  A(q)0f>*  and 

(Ml)  '  '  1^1*  <u  as**). 
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We  thus  arrive  at  the  final  form  of  the  approximating  system  for  (2.1)  in  XN(q)  as 
QNwN(t)  =  KNwN(t)+RNF(q,  t,  0NwN(t)),  I  >0, 

U‘i2>  wN(Q)^{QNy1RNu0(q). 

3.  Approximation  theorems  for  abstract  cystoma.  In  tins  section  we  shall  focus 
our  attention  on  the  condition  “q1  -*  4  implies  «<(/;£)”  of  Theorem  2.1(ii) 

and  present  results  on  the  convergence  of  solutions  of  the  approximating  systems  (2.7) 
to  solutions  of  (2.1).  We  state  and  prove  two  theorems;  the  first  is  applicable  to 
nonlinear  parameter  identification  problems  ($4)  while  the  second  will  be  used  in 
connection  with  linear  boundary  control  problems  in  }  5. 

Theorem  3.1.  Suppose  hypotheses  (H1HH3)  and  (H5)t(H9)  hold  and  let  qN,  q 
be  arbitrary  in  O  such  that  qA-*4-  Further,  suppose  that  |«o(<f * ) ~«o($)|^r •*  0  as 
N-*ao.  Then  the  mild  solutions  uN(t ;  qN)  of 

uN«)-Ar'(qN)uN(t)+F"(q",t,uN(t)), 

(3-1)  uN(0) = PN  (q^Uakq1*) 

converge  to  the  mild  solution  u(t;4)  of  (2.1)  for  each  fe[0,  T\  If  (t,v)-»F(4,t,v) 
is  continuous  on  [0,  T]*X,  then  the  convergence  lu^f;  qN)-u(f,  $)|,«  -»0  is  uniform 
in  t  on  [0,  T\. 

Proof.  Let  qN  -*  4  be  arbitrary  as  hypothesized.  Recalling  the  proof  of 
Theorem  2.1,  we  observe  that  one  easily  argues  existence  of  a  constant  K  such  that 
|uN(f;qv)|,"  S8T,  |u(f;4)|,"  SJC  for  all  AT  and  t  e  [0,  T],  Further,  we  see  that  for 
t  €  [0,  T\  we  have  (where  all  norms  are  H«") 

|nA'(l;«,',)-«(f;$)| 

S \T"(r,  <")P%W)M«K)  ~TN(f,  <j")P"(<j%o(4)l 
+  \TN(t ;  *")P*V)«o(4)  -  TN(f,  qN)uo{4)\ 
+\T"{r,q")uo(4)-nr,4)t*o(4)\ 

+  £ \TN(t-s,qf,)PN(qNmqN, s, u"(s;qs))-F(q*  s, u(s;4))} I ds 

+  £  |T*(r -s;  qN)PN(qNHF(qN,  s,u(si  4))~F(4,  s,  u(s;4)» I  ds 

+ £ \T"(t -s;  qN){PN{qN)-I}F(4,  *.u(s;  *))|  ds 

+£  \{TN(t-s-,q")-nt-s,4m4,  S,  UW;4))I  * 

-«t(N)+*2<A0+«s(A0 

4  £  |r"(t  s,  «"(*;  q"))-F(q",  s,  «(s;  4))W  ds 

+e4(A0+es(A0+ed(N). 

By  (Hi),  oar  hypotheses  mi  the  deftnitim  of  P"(q")  kt  (H3),  we  find  MWl* 
••  Alao,  |eKAOt*4Ke-TKPw(4>')-/)«o(4)l-*0  by 
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(H8)  and  (H7).  That  |«3(Ar)|  ■*  0  uniformly  in  /  on  [0,  T]  follows  directly  from  (H9). 
Moreover, 

(T 

\e4m\SMe“T  \F(qN,s,u(s-,d))-F(q,s,u(s-,q))\ds->0 
Jo 

by  (H6)(iv),  (H6)(iii)  and  dominated  convergence,  while 

|e5(JV)|S^e“'r  [  |{FV)-/}F(d,s,u(s;$))|ds-»0 
Jo 

by  (H7),  (H6)(iii)  and  dominated  convergence.  Finally, 

|e«(A0|  =  [ '  |{Tv(f  -  s ;  qN) -  T(t -s ;  )}F(*  s,  u  (s ;  4»|  ds  ->  0 
Jo 

by  (H9)  and  dominated  convergence  ((H8)  with  (H6)(iii))  for  each  fixed  t  €  [0,  T]. 

We  note  that  the  convergence  in  all  of  the  terms  above,  except  e6,  is  uniform  in 
t  on  [0,  T].  If,  in  addition,  the  continuity  hypothesis  of  the  theorem  obtains,  we  find 
that  {F(<J, 5,  u(s\ 4))|s e [0, 7"]}  is  a  compact  subset  of  X  and  the  convergence  in  the 
integrand  of  e6  is  uniform  in  t  and  s ;  hence  in  this  case  e6-»0  uniformly  in  t  also. 

We  have  thus  established  the  following  estimate: 

\uN(r,qN)-u(r,if)\ 

S  i  ei+Me’*  \'\F(qN,s,uN(s-,qN))-F(qN,s,u(s;il))\ds 
/-I  Jo 

*eN(t)+Me“T\  ki(s)\uN(s;qN)-u(s;<i)\ds 
Jo 

where  eN  -*0  as  N -*oo.  An  application  of  Gronwall’s  inequality  then  yields  that 
|u'v(f;qJV)-u(t;4f)|-*0  as  AT-*oo,  where  the  convergence  is  uniform  in  t  under  the 
added  continuity  hypothesis  of  the  theorem. 

Corollary  3.1.  Under  the  hypotheses  of  Theorem  3.1  ,uN  (t;q)-*u(,t\q)  for  each 
fixed  qeQ,  uniformly  in  t  on  [0,  T ]  if,  in  addition,  ( t ,  v)-*F(q,  t,  v)  is  continuous  on 
[0,  T]xAT. 

We  turn  next  to  convergence  results  needed  for  optimal  control  problems. 
Consider  for  fixed  q  e  Q  the  system 

ii(f)-A(<j)«(t)+G(q,f), 

(3-2) 

u  (0)  *  Mo, 

and  the  approximating  system 

uN{t) «  As(q)uf,(t)+Pf'(q)G(q,  t), 

«"(0)-F*(q),«o, 

where  G  has  the  form  0{a,  t)  -  y(q)cr(t).  We  assume  y (q)  e  T  where  T  is  a  subset  of 
<*-(*,•  with  X  a  given  subset  of  X.  We  further  assume 

<r  a  X,  X  a  given  subset  of  Lj(0,  T\ /?*). 

Theorem  3.2.  Assume  (H1MH3),  (H3),  (H7HH9).  Suppose  moreover  that  Jt 
is  compact  and  X  is  bounded.  Then  for  each  fixed  qeQ,  mild  solutions  us  of  (3.3) 
converts  to  the  mild  solutions  of  (3.2),  uniformity  in  cr  «  X,  y  «  T  and  t  a  [0,  T], 
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Proof.  We  consider  the  estimates  in  the  proof  of  Theorem  3.1  with  F(q,  s,  v)  = 
G{q,  s)  and  qN  =$=q  fixed.  Then 

k q ) - u  (t ;  q  )|  5  e2(Af) +  b3(N) + *,(JV) + e«(JV) 

(e  i  =  e« = 0),  where  e2  and  e3  are  as  before  and 

es=['\TN(t-s,q){PN(q)-imq,s)\ds, 

Jo 

86  “  1  *)l  <&> 

and,  as  usual,  all  X  norms  are  |*|,.  We  have  immediately  (using  (H7),  (H8),  (H9)) 
that  e2  and  e3  -►  0  as  Af  -»  oo,  uniformly  in  a,  y  and  /  e  [0,  T].  Also, 

\es(N)\  ZM.  e"T  fV" (q)-I)y(q)a(s)\  ds 
Jo 

%Mt-T  max  \(PN(q)-T)y,(q)\  (Vfc)!*. 

i  s/cm  Jo 

But  since  ^  is  compact  and  £  is  bounded,  this  latter  term  -+0  as  N  ■*  oo  uniformly  in 
•y  e r,  <r  e £  and  te[0,  T].  Finally, 

|es(N)|a  max  f  |{rN  (»-*;<*)- r(f-j;<t)}y.(q)IHs)|<6 

JO 

n.r 

){Tn  (/  -  j  ;  q )  -  T(t  -  s ;  q )}%  (q  )|2  di }  M 
o  J 

and  this  last  estimate  yields  e«(AT)  -*  0  uniformly  in  y  e  T,  <r  6  £  and  t  e  [0,  T],  again 
from  the  compactness  of  it,  (H9),  and  the  boundedness  of  £. 

As  we  have  previously  noted,  the  main  purpose  of  (H6)  is  to  allow  us  to 
guarantee  existence  of  solutions  of  (2.1)  and  (2.5)  on  fixed  finite  intervals  [0,  T] 
(see  Proposition  2.1).  The  condition  (H6)(iii)  is  used  in  the  proofs  of  Theorem  2.1 
and  3.1  only  to  establish  uniform  bounds  on  the  uN.  This  permits  us  to  employ  the 
local  Lipschitz  condition  (H6Mii)  and  to  appeal  to  the  dominated  convergence  theorem 
in  certain  arguments.  We  have  already  noted  that  (H6)(iii)  can  be  relaxed  to  “affine 
growth  at  oo”  (see  Remark  2.1).  With  an  alternative  approach,  one  can  relax  this 
growth  condition  even  further  and  still  obtain  die  conclusions  of  Theorems  2.1  and 
3.1  (with  the  other  hypotheses  remaining  unchanged).  Specifically,  for  N  sufficiently 
large,  the  initial  data  and  defining  operators  TN,  r  and  T  for  uN  and  u,  respectively, 
are  dose.  Thus  if  one  assumes  (in  place  of  (H6Xiii)) 

(A6)  (i)  For  each  qeO,  there  exists  a  solution  u(t;  q)  of  (2.1)  on  [0,  T\  and 

(ii)  There  exists  *3eZ,2(0, T)  such  that  |F(fl, t, 0)|« Sk3(t),  for  q,qeQ  and 

fe[0,n 

it  is  rather  tedious  but  not  difficult  to  show  that  for  N  sufficiently  large,  all  uN  defined 
by  (2.7)  exist  on  [0,  T]  and  lie  in  some  bounded  neighborhood  of  u,  the  solution  of 
(2.1).  (The  arguments  involve  use  of  classical  fixed  point  ideas  to  obtain  solutions 
on  some  interval  [0,Ai]  and  then  continuation  to  [6i,2dt],  •  •  •  ,etc.)  The 
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condition  (A6)(ii)  can  then  be  combined  with  (H6)(ii)  to  obtain  domination  of 
terms  such  as  F(q,  s,  uN(s\  q)).  Thus  ail  of  the  arguments  behind  Theorems  2.1 
and  3.1  remain  valid,  the  hypotheses  being  changed  only  in  that  (A6)(i),  (A6)(ii) 
replace  (H6)(iii),  and  the  conclusions  changed  only  in  that  they  can  be  obtained  only 
for  N  sufficiently  large  (which  of  course  is  not  important  theoretically  in  approximation 
results  such  as  those  discussed  in  this  paper). 

Regarding  the  assumption  (A6)(i),  we  note  that  there  are  various  conditions  one 
might  impose  on  F  to  insure  existence.  For  example,  monotonicity  hypotheses  might 
be  assumed  so  that  -(A  +F)  is  maximal  monotone  and  one  could  then  appeal  to 
standard  existence  results  [9],  [19].  In  §6  we  present  numerical  results  for  our 
approximation  scheme  for  an  example  (Example  6.5)  in  which  (H6)(iii)  is  not  satisfied, 
yet  (A6)(i)  and  (A6)(ii)  do  hold.  However,  we  shall  not  pursue  any  of  the  theoretical 
ideas  here,  since  this  is  really  not  the  focus  of  our  presentation. 

4.  Examples:  Parameter  Identificatiou  fas  hyperbolic  and  parabolic  equations. 

We  turn  now  to  an  application  of  the  results  developed  in  the  preceding  sections  to 
identification  problems  for  specific  equations.  A  fundamental  requirement  in  both 
Theorems  3.1  and  3.2  is  that  the  conditions  of  (H7),  (H8)  and  (H9)  be  verified.  As 
we  have  indicated  earlier,  the  convergence  statement  of  (H9)  can  be  obtained  rather 
easily  for  our  schemes  from  (H7),  (H8)  and  some  standard  approximation  results  from 
linear  semigroup  theory.  We  state  here,  for  our  future  reference,  one  version  (due  to 
Kurtz  [24])  of  these  approximation  theorems. 

Proposition  4.1.  Let  (91,  |  •  |)  and  (9IN,  |  •  |w),  N  =  1, 2,  •  *  • ,  be  Banach  spaces 
and  letirN  :9-*9N  be  bounded  linear  operators.  Assume  further  that  3~(t)  and  it) 
are  linear  Co-semigroups  on  91  and  9tN  with  infinitesimal  generators  si  and  siN, 
respectively.  If 

(i)  lhri*-.oo  |wnz|n  « |z|  for  all  z  €  91, 

(ii)  there  exist  constants  M,  £  independent  of  N  such  that  \&N (t )|  SM<",  f  £0, 

(iii)  there  exists  a  set  9^9  such  that  2 <=  Dom  (.si),  and  (Aa-jtj£~9 

for  some  ho>0, 

(iv)  for  z  e9  we  have  hm*... \sisirNx  - wNstz \s  - 0, 

then  lim  1 3TN(t)srNz  -irA,5'(f)zU  -0  for  z  6  91,  uniformly  in  t  on  compact  subsets  of 

[0.®). 

We  note  that  the  requirement  in  (iii)  implies  that  9  is  a  core  [21,  p.  166]  of  st\ 
this  is  easily  seen  using  the  fact  that  si,  being  an  infinitesimal  generator,  is  dosed  and 
ifj -si)'1  k  bounded  for  A  sufficiently  large.  The  proposition  then  follows  directly 
from  Theorem  2. 1  of  [24]  taken  with  subsequent  remarks  [24,  p.  361]  of  that  reference. 
Obviously,  (iii)  in  our  statement  above  could  be  replaced  by  die  hypothesis  that  9 
be  a  cose  of  si.  Further,  we  remark  that  the  requirement  fS  *9  is  superfluous  in  (iii) 
if  one  verifies  that  (A0-  si)&  »  91  for  Ao  *p(si).  In  this  case  one  can  easily  demonstrate 
directly  that  9  is  a  core  tot  the  generator  si. 

In  the  examples  we  discuss  below,  we  shall  use  the  notation  A,  F,  AN,  P",  FN 
to  denote  die  specific  operators  in  each  example,  since  this  will  facilitate  reference 
back  to  the  basic  theorems  of  If  2  and  3  and  should  cause  no  confusion  for  readers. 
However,  we  shall  adopt  distinct  notation  for  the  various  state  spaces  X,  XN  within 
the  context  of  each  example. 

Example  4.1.  Hyperbolic  equations.  We  consider  the  one -dimensional  hyperbolic 
equation 


(4.1) 


e»  +fi*o,+fi»v  +/(«*,  t,  x,  c,  v,) 
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with  initial  and  boundary  conditions 


(ICO 


r(0,x)  =  I  q*<h(x), 

i-i 

t>,(0,jr)=  £  q‘s4> f(x)  forOSxSl, 

i-i 


(BCO  v(t,0)  =  v(t,  1)  =  0  for  f  >  0* 

where  c  qteRm  and  the  remaining  q,,q\  are  scalars.  The  vector  parameter 

q  of  i  2  thus  has  dimension  *  *  3m  +3.  (If  the  output  maps  such  as  V  or  C  in  the 
fit-to-data  functions  (2.3)  and  (2.4)  depend  explicitly  on  some  parameters  q2 ,  •  •  • ,  q?, 
we  assume  with  no  loss  of  generality  that  these  have  been  embedded  in  the  (or  q« 
or  qs)  vector.) 

We  remark  that  we  do  not  formulate  nontrivial  boundary  conditions,  possibly 
depending  on  parameters,  in  (4.1HICiHBCi);  however,  it  is  easily  seen  that  by 
simple  transformations  such  generalities  actually  can  be  included  in  our  formulation 
above.  Consider,  for  example, 

(4.2)  ott-qxv„ 

with  initial  and  boundary  conditions 

v(0,x)~q4+(x), 

o,(P i,x)-qstt>(x), 

v  (r,  0)  -  <f  761(f),  v(t,  1)»  qj>i(t), 

where  61,  b2  are  twice  continuously  differentiable  functions.  Employing  the  standard 
transformation  w(r,j0  =  »(r,x)-(l-x)q76i(f)-xqs6a(f)»  we  find  that  (4.2HIC2)- 
(BCj)  can  be  reformulated  as  a  special  case  of  (4.1HICxMBCi). 

Returning  to  (4.1),  we  proceed  in  the  usual  manner  to  rewrite  the  equation  with 
boundary  and  initial  conditions  as  an  abstract  evolution  equation.  Let  A  denote  the 
Lapladan  operator  l?/dx2  in  H° - L2( 0, 1;  R)\  here  and  below  the  Sobolev  spaces  H1 
consist  of  R1 -valued  functions  on  [0, 1]  taken  with  their  usual  inner  product;  unless 
otherwise  specified.  It  is  well  known  that  A  with  Dom  (A)-Hj  (Iff2  is  a  self-adjoint 
operator  in  Hc  satisfying  <-Ar,  r)s|z|2  for  all  z  e  Dom  (A).  We  impose  the  following 
additional  assumption  on  the  coefficient  qx  in  (4.1): 

(HQ)  There  exist  positive  numbers  q  \  and  q  \  such  that  q  e  Q<=Rk  implies  q  fSqiS 

qY. 

For  a  given  q  e  Q<=Rim*3  we  of  course  mean  by  q  1  the  first  coordinate  of  the  vector 
«-(«!.•••.  qt)  where  q,  -  (q},  •■ ,  q?),  i  »  4, 5, 6. 

Having  thus  assumed  hypothesis  (HQ)  for  a  given  fixed  parameter  restraint  set 
Q,  we  endow  the  set  Ho  with  a  family  of  inner  products 

<m  «>.-£  4.  £  £-<**„..> 

where  q  ranges  over  Q.  In  view  of  (HQ),  the  space  (Ho,(*,  •)„)  is,  for  each  qeQ,  a 
Hilbert  space  which  we  denote  by  V(q0.  The  space  X(q)  of  f  2  is  chosen  for  this 
example  to  be  Xiq)-X(q)»V(qi)xlT  with  the  usual  product  topology  generated 
by  <(wi,  wj),  (z u  **)>«  -  <wi,  ri)fl  +<wa,  *j).  Condition  (H2)  is  obviously  satisfied  since 


(IC2) 

(BCj) 
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for  any  q,qeQ  we  find  \z\fSX\z\q  for  all  z  eX  where  X~(qi/qi)i/2.  Introducing 
the  variable  w(i)~vB  we  may  rewrite  (4.1HICiHBCi)  in  X  by 


(4.3) 


t>0. 


\w(0)/  l£  fW' 


where  (4«,  0() 6  Dom  (A(q))  =  Ho  0 H1  x f/J, 

0 


and 


F(q,f,t>(r),  w(f))“(,,  ,  /?  v  UY 


Before  turning  to  a  careful  discussion  of  (4.3),  we  define  the  operators,  etc.,  needed 
in  formulating  the  modal  approximation  scheme  associated  with  (4.1). 

Since  the  operator  -A  is  self-adjoint  with  compact  resolvent;  standard  results  in 
spectral  theory  and  the  theory  of  bilinear  forms  (see  [15,  p.  1331],  [39,  p.343], 
[33,  pp.  247-254])  are  applicable.  Defining  «ty(jr)  ■  (V2//ir)  sin jwx  and  ♦/(*)» 
•Jl  sin  jwx,  we  find  that  {£;}“  i  and  {♦/}“  i  constitute  complete  orthononnal  sets 
(CONS)  for  V(l)(«Ho)  and  H°,  respectively.  The  corresponding  modal  subspaces 
Xs  (q)  -X*(q)  of  X(q)  are  then  defined  by  Xs (q)  -  span  { ’fit, ,  •  •  • ,  0&,)  where 

It  is  easily  seen  that  U  forms  a  CONS  for  X(q*)  where  q*  -  (1, 0,  •  •  • ,  0), 

and  a  complete  orthogonal,  but  not  normal,  set  for  X (q),  for  q  « (q»,  •  •  • ,  q*)  with 
qi>0,  qi  *  1.  We  note  also  that  $<,  are  eigenfunctions  of  A  corresponding  to  the 
eigenvalues  A  1, 2,  ■  •  • 

The  subspaces  X"(q)  and  the  corresponding  orthogonal  projections  PN(q)  (see 
(H5))  having  been  thus  defined,  the  modal  approximation  operators  Aw(q)  for  A(q) 
are  determined  as  in  (H3).  The  corresponding  matrices  QN  and  KN  of  (2.10),  which 
in  this  case  are  x  2N-  matrices,  are  readily  seen  to  be  given  by 


(4.4)  Qn  ■  (Hag  (qlt  •  •  • ,  qi,  1,  •  •  • ,  1), 

where  the  qt  and  1  each  appear  N  times,  and 


where  the  D*,  /  —  1, 2, 3,  are  NxN  diagonal  matrices  defined  by  D" >diag(wqi, 
2wqi,  •  •  • ,  Nvqi),  D"  •  diag  (q»/w,  q»/2ir,  •  •  • ,  qt/Nv) -Du  and  D*  •  qj.  Recall¬ 
ing  (2.9),  we  observe  that  in  this  case  the  projection  operators  F*(q)  are  actually 
Independent  of  q. 

We  are  bow  in  a  position  to  verify  that  (HI)  and  (H9)  obtain  for  die  hyperbolic 

TKBOUM4.1.  Assume  that  (HQ)  holds  and  kt  (qw)  be  an  arbitrary  sequence  in 
QcJtim*s  converting  to  («<}.  Then  the  operators  A(f)  and  A*(q*)  of  (4.3)  and 
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the  associated  modal  approximations  described  above  generate  Co-semigroups  T(t;  4) 
and  on  X(f)  and  \Ts{t-,qN)z~T{t\4)zU”-*0  for  z  eX(4),  uniformly  in 

t  on  compact  subsets  of  [0,  oo).  If  we  further  assume  (H4)  (Q  is  compact),  then  there 
exists  tteR,  independent  of  q  eO  and  N,  suck  that  \T(t\q)\Se"^  and  |TM(r;4)|Se"< 
for  all  qeQ,  tiO,N  •  1,2,  •  •  • 

Proof.  For  any  q  e  Q,  a  straightforward  calculation  shows  that 

“•W’Ci  o) 

with  T)om{A<Jlq))~H\  Dfi'xHo  is  a  symmetric  operator  in  X(q).  Furthermore, 
Dom  (A0(q))  is  dense  and  Ao  is  invertible.  It  follows  [28,  p.  97]  that  A0  is  skew  adjoint 
and  hence  by  Stone’s  theorem  [26,  p.  252],  [41,  p.  345]  generates  a  unitary  group  on 
Defining  the  operator  A(q)  on  X(q)  by  A(q)(zu  zi)  -  (0,  q$z i +qiZi),  it  is  easily 
seen  that  A(q)  is  bounded.  Indeed,  using  the  fact  that  the  Ho  norm  is  stronger  than 
the  H°  norm,  one  finds  \A(q){zu  z2)|,  Sc(qu  qt,  qj)|(zi,  z2)|,  where  die  constant  c  is 
bounded  above  uniformly  on  Q  if  (HQ)  and  (H4)  hold,  say  c(qu  q2,  qj)So>. 

We  thus  find  that  A(q)  ~Ao(q)+A(q)  is  the  perturbation  of  A0  by  a  bounded 
operator  and  hence  [21,  p.  495],  [30,  p.  80]  generates  a  CD- semigroup  T(r;  q)  on  X{q) 
satisfying 

(4.6)  |jf*(t;4)|aexp{c(g1,g2,gs)r], 

where  there  exist  o>  >0  independent  of  q  such  that  c(q i,  q2,  qs)  Set  in  case  (HQ)  and 
(H4)  obtain  (or  in  case  (HQ)  holds  and  q  lies  in  a  bounded  subset  of  Q). 

As  we  have  pointed  out  earlier  (see  die  remarks  in  i  2),  AN(q)  is  a  bounded 
linear  operator  for  each  N  and  hence  generates  a  Co-semigroup  on  X(q).  Assuming 
that  (HQ)  holds  and  q  lies  in  a  bounded  subset  of  Q,  we  have  that  A(q)  is  the 
infinitesimal  generator  of  a  Co-semigroup  satisfying  |T(r;g)|SeM  so  that  A(q)~toI 
is  the  generator  of  a  contraction  semigroup  and  is  hence  dissipative  [22,  p.  90],  [30, 
pp.  14-17].  That  is. 


(A(g)*,*>  Sv(:,r) 

for  all  z  e Dom  (A(q)).  From  the  definition  of  Ax  it  follows  that  for  z  e  X(q) 

(4.7)  (Aw(g)z,z), -(Alq^z,  P*z >,  S«|J>wz|;s«l*lJ. 

since /’"(qlisthe  orthogonal  projection  of  JF(g)  onto  #*({).  Hence  we  findlT^fo^ls 
e*,  as  desired  in  the  second  condushm  of  the  Qieorem. 

We  make  use  of  Proposition  4.1  to  establish  the  convergence  results  of  the 
theorem.  We  take  for  our  discussions  &-X{g)  and  St N -X(q"),  which  of  course 
satisfy  die  conditions  of  (H2).  The  above  arguments  yield  immediately  that  (ii)  of 
Proposition  4.1  holds  for  our  family  of  semigroups  TN (t;  q").  Letting  w*  be  the 
canonical  isomorphism  from  X(<f)  to  X(qN),  we  see  immediately  from  the  hypothesis 
qM  ■*  4  and  the  definition  of  the  norms  In  X{q)  that  |w*z  (,» |z  |*  so  that  (i)  is  satisfied. 

We  define  and  have  at  once  that  $eDom(A(4))  and  J- 

X(4).  From  the  definition  of  A(f),  the  tact  that  Jtl-A(4)  b  invertible  tor  A 
sufficiently  large  and  that  the  4>/  are  eigenfunction  of  A,  it  is  easily  argued  that 
(A/ -A(4)ya>  -  9  so  that  {k  -Amb-Xtf)-,  hence  (iii)  is  satisfied. 
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Finally,  suppressing  the  notation  it*  *  for  the  canonical  isomorphism  (see  our 
comments  in  S  2),  we  have  for  each  z  e  ®  and  AT  -  Af(z )  sufficiently  large  (then  PNz  =  z ) 

|A  V)z  ~A(4)z\ « \PNA(qN)PNz -A(4)z\ 

SlJ»*,A(flw)z  -P^AdW+^AWz  -A(j)z\ 

S  |A(q")z  -i4(4)z|+|(P* -I)A(4)z\ 

where  all  norms  are  j  ■  I,*.  Since  qN  -*H  by  hypothesis,  the  form  of  A(q)  yields  that 
the  first  term-»0  as  N-*<x>.  From  the  completeness  of  the  {0*}  in  X(q*)  (see  our 
remarks  above),  the  equivalence  of  norms  and  thus  the  strong  convergence  PN  ■*  I  in 
any  of  the  norms,  we  obtain  that  the  second  term-»0.  This  completes  the  proof  of 
Theorem  4.1. 

We  return  to  the  abstract  nonlinear  equation  (4.3)  to  consider  conditions  on  /  of 
(4.1)  so  that  F  of  (4.3)  will  satisfy  (H6)  and  hence  the  results  of  $S  2  and  3  will  be 
applicable.  Define  <?*  -  {q«  e  Rm |q  €  0}  where  Q  cR 3m *3  is  given  for  (4. 1).  We  impose 
the  following  hypotheses  on  /. 

(H6*)  The  nonlinear  function  /:  06 x [0, T]x[0,  l]x/f  xR  -*R  satisfies: 

(i)  For  each  (qs.u,  w)eQ*xHJ,  the  map  (t,x)-»f(q6,t,x,  v,w)  is 
measurable. 

(ii)  For  each  constant  M >0,  there  exists  a  function  fci  *  ki(M)  in  L2(0,  T) 
such  that  for  all  (q*.  t,  x)e  Q*x[0,  T]x[0, 1]  we  have 

1/(4*.  t,x,vu  w)-f(q6,  t,  x,  v2,  w)|  S £i(f)|ui  -  u2| 

for  all  (vi,  w)eK3  with  \v,\sM. 

(iii)  There  exists  a  function  k2  in  L2(0,  T)  such  that 

1/(4*.  t,  x,  v,  H>t)  -/(q*,  f,  X,  v,  Wi)|  S  £i(t)| Wj  -  w2| 

for  all  (q*.  t,x)e  O«x[0,  T]x[0, 1],  and  v,  Wi,  wteR. 

(iv)  There  exists  a  function  £»  in  L2(0,  T)  such  that 

1/(4*.  i,  x, ».  o)|  a  £j(fM|»| + 1} 
for  aU  (q*.  t,  x,  »)«  Q*x  [0, 7}x  [0,  l]xlf. 

(v)  For  each  (t,x,  v,  w)€[0,  T]x[0,  l]x/fJ,  the  map  q*-*/(q*,f,x,u,  w)  it 
continuous. 

Employing  rather  standard  arguments  and  results  from  analysis  (e.g„  see  [15, 
Lem.  16(b),  p.  196]  in  connection  with  (i)),  it  is  quite  straightforward  to  verify  under 
(H2)  that  (H6*)  for  /  implies  (H6)  for  F  in  the  example  under  consideration.  We  can 
therefore  appeal  to  Proposition  2.1  to  guarantee  existence  of  a  unique  mild  solution 
of  (4  J)  for  any  q  c  Q. 

Summarising,  we  have  shown  that  under  (HQ),  (H4)  and  (H6*)  the  conditions 
(H1MH9)  hold  for  the  abstract  formulation  (4.3)  of  (4.1HIC1HBCi)  when  consider, 
ing  the  modal  approximation  scheme 


(4.8) 


sC>)) ‘  '“«»•  *"«'»• 


in  JF*(q).  The  convergence  of  Theorems  3.1  and  3.2  is  thus  assured  and  we  may, 
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when  an  appropriate  fit-to-data  function  is  defined,  appeal  to  Theorem  2.1  to  obtain 
a  solution  o f  the  associated  identification  problem  for  (4.IHIC1HBC1).  For  example, 
suppose  we  me  given  observations  ${U) €*',  i  - 1,  •  •  • ,  r,  for  (0(4,  x»),  •  •  • ,  v(lh  xt)) 
in  (4.1)  where  05fi<  •  •  •  <trST.  Let  u(f;4)«(ui(r, •;«).  «a(t.*;q))  denote  die 
unique  mild  solution  of  (4.3)  where  we  observe  that  ui(t ;  q)eHl  for  each  <,4,  so  that 
potntwise  evaluation  in  [0, 1]  is  a  meaningful  operation.  Let )  have  die  form  given 
in  (2.4)  where  now  4(»bq)-col(u1(fbx1;q),---,«i(f,,Xi;q))  and  C(thq)is  an  /x 
1-  matrix  depending  continuously  on  q  in  Q  for  each  i.  Hum  dearly  this  /  satisfies 
the  continuity  requirements  of  Theorem  2.1.  Furthermore,  H  is  easily  seen  that  die 
initial  data  in  (4.3)  depend  continuously  on  q.  For  the  modal  approximations,  recall 
that  P*(q)  is  independent  of  q  and  finally  note  that  the  continuity  of  q •*  Tfl(t\q)z 
follows  directly  from  the  forma  of  As(q),  QN,  and  KN  given  explicitly  in  (2.10),  (4.4) 
and  (4.S).  The  conclusions  of  our  deliberations  for  Example  4.1  may  thus  be  dated: 

Theorem  4.2.  If  (HQ),  (H4)  and  (H6*)  obtain,  then  the  problem  (1DN)  for  (4.8) 
with  J  as  defined  above  has,  for  each  N,  a  solution  4*  e  O  cR3m*3.  Letting  (4*'}  be 
any  subsequence  of  (4*}  converging  to  ijeQ,  then  4  is  a  solution  of  the  problem  (ID) 
for  (4.3)  and  moreover  for  each  t  e[0,  T], 

|(t>w'(t;  4*0,  w\r,4N‘))-(ut(f,4),  «a(r;4))l-0 

as  Nj-*<x>  where  («i,  u2)  is  the  solution  of  (4.3)  and  the  norm  is  that  of  Hi  xH°. 

Remark  4.1.  We  remark  that  the  dependence  of  the  norm  on  q  in  the  above 
treatment  of  hyperbolic  systems  is  somewhat  artificial.  While  one  cannot  effectively 
rescale  the  time  variable  to  remove  the  qt- dependence  in  problems  where  sampling 
times  (observations)  are  important,  one  can  rescale  die  state  variables  (use  w(r)« 
l/Vqiet  in  place  of  die  variable  used  in  (4.3))  til  avoid  use  of  a  weighted  norm.  We 
chose  not  to  do  that  in  our  computations  for  Example  4.1.  A  preliminary  consideration 
lends  one  to  conjecture  that  such  a  rescaling  does  not  result  in  simplification  from  a 
numerical  viewpoint. 

Example  4.2.  FaraboUc  equations  I.  For  our  second  class  of  examples  we  discuss 
scalar  parabolic  equations 

(4.9)  e,  - ^ (p»,),  +q2v  +f[q\,  •  •  • ,q?,t ,  x,  v) 

for  f  >  0,  x  e  [0, 1]  with  initial  and  boundary  conditions 

(IQ  e(0,x)»  £ 4sdi(*).  OSxSl, 

1*1 

(BQ  */»«,•) -0,  /■  1. 2. 

Here  we  assume  that  ^eH0,  »(r,x)  (or  u(»,x;q))  is  in*,  and  q  •  (41,42. 4s.  4s)  <Q 
where  <?cgJ*+1  and  q,-(q/\  •  •  •  ,4")  for /-3,4.  The  operators  *i,  *j  defining 
the  boundary  conditioas  have  domain //2  sod  are  given  for  by 

(4.10)  Rfb -allib(Q)+ait+'(Q)+alit{l)+aiob,{,l), 

where  a»u*.  We  impose  die  fotiowing  conditions  on  k  and  pin  (4.9)  and  the  a*: 

(HP1)  The  functions  p.p,  and*  are  continuous  with  k(t)>0,p(x)>0for0SxSl. 
(HP2)  The  matrix 

Cu  #i>  «u\ 

u  *»  «»  ««/ 

hat  rank  2  and  wehavepfOHanaja-aixosiJ-pdHaijOM-awaas}. 
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We  shall  also  assume  that  our  parameter  set  0  satisfies  hypothesis  (HQ)  already  given 
above. 

To  apply  the  results  of  S8  2  and  3  we  again  rewrite  our  system  as  an  abstract 
evolution  equation.  Although  we  consider  here  only  a  scalar  equation,  the  results  we 
present  generalize  immediately  to  coupled  systems  of  parabolic  equations.  Indeed, 
we  present  numerical  results  in  i  6  for  examples  of  such  systems  which  are  of  interest 
in  biological  applications  (see  [4]),  among  others. 

We  define  the  general  Sturm-Liouville  operator  A(q)  (our  operator  here  is  the 
negative  of  die  one  found  in  the  literature  cited  below)  in  H°  by  Dom  (A(q))~ 
{#  e H2\Rfi>  -0,  /  - 1, 2}  and  A(q)#  -  *'l(q,p*A  +q2*  Then  (4.9MICMBC)  can 
be  rewritten  as 

v(t)  * A(q)vU)+F(q ,  t,  c(f)),  t>0, 

®(0)»  £  h 
i-t 


where  F(q,  r,u(r))=/(q4,r, •))  and  the  equation  is  taken  in  the  state  space 
2f(q)*Sfr»(ff°,<-,  •))  with  <4J,«fr)*J l  d(x)d(x)k(x)  dx.  We  note  that,  unlike  in 
Example  4.1,  here  the  state  space  actually  doesn’t  depend  on  q. 

Spectral  results  for  the  operator  A(q)  are  readily  found  in  the  literature — e.g., 
see  [28,  p.  182],  [20,  p.  126].  The  hypotheses  (HP1),  (HP2)  imply  that  A(q)  is  self- 
adjoint  and  its  spectrum  consists  of  a  countable  number  of  real  eigenvalues  (A,(q)}/“  i, 
each  of  multiplicity  less  than  or  equal  to  2,  and,  moreover,  these  eigenvalues  can  be 
ordered  so  that  -co<  •  •  •  s»A,5A,_iS  •  •  •  SAt<oo.  The  eigenfunctions 
corresponding  to  A(q*)  where  q*»(l,  0,  •  •  • ,  0)  form  a  CONS  in  £. 

We  further  observe  that  the  eigenvalues  (A,(q)}  of  A(q)  are  bounded  above 
uniformly  in  q  on  bounded  subsets  of  Q.  This  is  easily  seen  as  follows:  Let  X*  be  the 
eigenvalues  of  A0**  A(q*)  (i.e.,  Aod-A-1^,),)  with  corresponding  eigenfunctions 

From  our  remarks  above  we  have  X; <#,/ *  1,2,  •  •  • ,  for  some  positive  finite 
constant  <3.  The  eigenvalues  for  A(q)  are  then  found  to  be  A,(q)«>qiXj+q2  (with 
eigenfunction  %)  so  that  we  find  A/(q)S«  on  bounded  subsets  of  O,  where  a>  is 
independent  of  q  (but  depends,  of  course,  on  the  particular  bounded  subset  of  Q 
involved). 

We  define  the  approximating  modal  subapaces  of  X(q)«$  by 
span  {'Pi,  •  •  • ,  ¥*}  and  let  P" :  JP*  denote  the  associated  canonical  orthogonal 

projections.  As  before,  determine  the  operator  AN(q)  and  F*  by  A(q)« 
PsA(q)Ps  and  F*  -P*F.  We  have  die  following  convergence  results. 

Thbomm4J.  Suppose  that  (HQ),  (HP1)  and  (HP2)  hold  and  qff,4eQcRim*1 
an  sack  that  qH  -*4  at  N-+  ao.  Then  A(4)  and  A*(q ")  generate  Co-semigroups 
T(f;  4)  and  T"(f\  f  )  on  it  and  |r*V}  q "jg  -  T(t\ 4)t\*0fort  e  it  with  the  conver¬ 
gence  uniform  in  t  m  compact  mbeets  of  [0, »).  Furthermore,  4  (H4)  obtains,  then 
then  exietto  conetontmindepondontofNond  qmeeh  duu  |T(l;  #)j  and  |rN(f;  q)\  * 

o'"  for  t  >0, 4 a 0,  and  N  ■  1, 2,  •  •  • 

Proof.  Let  O  be  any  bounded  subset  of  Q.  Then  our  remarks  above  imply 
i»tstsaniafJ-w«?)sach  that  the  sod  edjoint  opereaorA(q),qeC?  has  its  spectrum 
lying  in  (-<©,#).  Hence  (see  [36,  p.  349])  A(q)-AI  h  dlwlpetive,  i.e.,  <(A(q)- 
«!/)*,  r>*0  for  al  e*Dom(A(q))  and  q  a  4.  For  A0«r(A(q)),  we  have  [28,  p.  180] 
that  A(q)—Af  km  compact  lesoheat  so  that  In  particular  wa  have  (A(q)- 
AI)  Dorn  (A (q)) -  #  for  A  >0  property  cboaen.  It  foBowt  Immediately  [30,  p.  17],  [2, 
p.  173],  [22,  p.  87]  that  A  (q)  -Jilt  maximal  dhslpatire  snd  gaimstes  a  Cp-seatigroup 
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of  contraction*  on  ft.  Arguing  as  in  (4.7),  we  have  that  As(g)-tSI  is  dissipative, 
uniform  in  q  €  <5  and  N  *  1, 2,  •  •  •  (it  is  maximal,  being  defined  on  all  erf  ft — see  [22, 
p.  86]),  and  hence  is  also  the  generator  of  a  Co-semigroup  of  contractions. 

The  above  remarks  obviously  apply  if  we  choose  0  =  where  qN  -*q 

or  if  0  m  Q  where  Q  itself  is  compact,  i.e.,  (H4)  holds.  To  obtain  the  convergence 
results  of  the  theorem,  trim  (3  ■*{qJlU(<f}  and  use  Proposition  4.1.  Here  =  # 

and  coaditioRs  (i)  and  (ii)  of  the  proposition  dearly  are  satisfied.  Let  S'  =  US-i  ftN 
so  that  §cDom  (A(4))  and  9  is  dense  in  ft.  From  our  remarks  above,  we  have 
A^)*/- Ay ($)♦,.  For  A>max;A/(4)  we  have  A I-A(4)  invertible  where  (A I- 
A(4))*,  -  (A  -  A /(4))%  mi  it  follows  that  (A/  -  A(4))$  =  9.  Hence  (iii)  of  Proposi¬ 
tion  4.1  is  satisfied.  The  arguments  that  (iv)  is  satisfied  are  exactly  analogous  to  those 
used  to  complete  the  proof  of  Theorem  4.1,  here  the  completeness  of  the  {¥,}  yielding 
PN  -*I  strongly  in  ft.  We  thus  have  established  the  results  of  Theorem  4.3. 

We  turn  finally  to  conditions  on  /  in  (4.9)  that  will  insure  that  F  of  (4.11)  satisfies 
(H6).  Let  Qtm{q*eRm\q  e  <?}  where  Q  is  a  given  subset  of  R2m+i. 

(H6**)  The  nonlinear  function /:  Q4  x [0,  T]x[0,  l]xR  -*R  satisfies: 

(i)  For  each  l<j4,  v)eQAxR,  the  map  ( t ,  x)-*f(q *,  t,x,  v)  is  measurable. 

(ii)  There  exists  a  function  £»  in  L2(0,  T)  such  that  |/(q4, 
f(q*t  t,  x,  c2)|  S  ifi(r)|oi  -  e2|  f«  all  (q*.  t,  x,  v,)  e  Q*  x  [0,  T]  x  [0, 1]  x  R. 

(iii)  There  exists  a  function  Jt2  in  L2([0,  T]x[0,  ID  such  that  |/(q«,  t,  x,  0)|S 
i2(t,  x)  for  all  q4  e  0«- 

(iv)  For  eadi  (r,  x,  v)  in  [0,  r]x[0,  l]xH,  the  map  q4-»/(q4,  t,  x,  u)  is  con¬ 
tinuous  on  Q4. 

It  is  an  easy  exercise  to  verify  that  (H6**)  for  f  implies  (H6)  for  F  (note  that  in 
this  example  die  condition  (H2)  is  trivially  satisfied).  We  thus  may  invoke  Theorem 
3.1  for  convergence  of  our  modal  approximation*  defined  in  by  the  equation 

dw(t)  -  A"(q)v  N(t)+PNF(q,  t,  vN(t)), 
tAO)-/*  1  q%. 


For  dieee  parabolic  systems  defined  in  H°,  the  question  of  an  appropriate  cost 
functional  is  somewhat  mote  delicate,  since  in  general,  point  evaluation  may  not  be 
meaningful  Hue  pnmltiUtly  fa  iUfftirunt  arrrnnrh  trill  Its  itisnmrril  brl — 1  ir  r  -hirer 
a  coat  functional  /  as  in  (2.3)  where  now  «  -  v  is  the  mild  solution  of  (4.11)  end  we 
might  assume,  for  example,  that  (x,q)-*  Y(tkx,q)  is  continuous  for  each  rt.  We  are 
again  in  a  position  to  employ  Theorem  2.1  (taking  info  account  the  comments  in 
Remark  2.2)  to  establish  the  foOowtng  results  for  this  example. 

Thboubm  4.4.  Suppose  (HQ),  (H4),  CH6**)  and  (HP1),  (HP2)  an  satisfied.  Then 
the  problem  (ZD*)  for  (4.12)  with  J as  In  (2.3)  has  a  solution  a  <?  c  Ra"*J  for  sack 
AT  - 1, 2,  •  •  • .  (t"1)  Is  any  subsequence  of  {4")  converging  toiuQ,  then  4  is  a 

solution  of  (ID)  for  (4.11)  and  for  sack  t  ntO,  T]  we  have  l»'(f;4N')-c(f;4)|-*0,  as 
N,-*ao,  when  v,  v*1  an  die  mild  solutions  of  (4.11),  (4.12)  rtsptctioely  and  die  norm 
Is  that  of  N*. 

Example  4  J.  Farabottc  equations  D.  We  consider  again  dm  parabolic  equation 
(4.9)  with  bridal  condition  (1C)  and  boundary  condition  (BC)  but  with  ntiffatly  mote 
restrictive  conditions  on  tha  boundary  operators  R  than  those  given  in  (4.10)-(HP2). 
Wa  treat  problems  with  boundary  opera  ton  chosen  ttum  the  standard  term- 
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Lkraviflc  operators  (see  [14,  p.  143],  [13,  p.  291}) 

(A)  *,*«*<))  or  *^=#(1); 

<B)  or  *^-*'(1); 

'  (C)  tft*»#'(0)-<ri*(0)  or  /lj#-^(l)+<r^(l),  <rt,<r,>0; 

(D)  /?,^  =  ^(0)-4r(l)  and  *,d-p<0)#'«»-p<l)d’(l). 

In  the  sequel,  in  referring  to  the  standard  Sturm-LiouviUe  boundary  conditions 
(SLBC),  we  shall  mean  any  combination  of  chokes  for  Ri  and  Ra  from  (4.13)  (AHC) 
or  the  periodic  boundary  conditions  arising  from  the  choice  of  Ru  R2  given  in  (4.13) 

(D). 

Our  main  objective  is  to  discuss  the  use  of  the  simple  pointwise  fit-to-data  criteria 
(4.14)  Ji(<M>(-;4).y)-  t  im-CfoqMM)? 

i-t 

as  defined  in  S  2  (see  (2.4)  and  the  discussions  thereof).  When  treating  (4.9)  in  £  as 
we  did  in  Exan  pie  4.2,  it  is  by  no  means  dear  that  the  associated  ID  problem  with 

(4.14)  is  well  posed.  Indeed,  one  must  first  justify  die  pointwise  evaluation  (in  the 
spatial  coordinate)  of  e  involved  in  defining  f ;  assuming  that  this  can  be  done,  one 
must  entertain  a  second  difficulty  in  that  die  convergence  (of  Theorems  3.1  and  4.3) 
of  vs{f,qs)  to  t>(f;4)  is  in  the  3?  (i.e.  H°)  norm.  Since  /t(*,  •,  f)  is  not  continuous 
on  QxC(Q,  T;  jt),  the  results  of  Theorem  2.1  are  not  directly  applicable. 

We  turn  (fast  to  the  difficulty  raised  by  point  evaluations  in  (4.14).  hi  this  regard 
we  note  that  the  mapping  w-+Ji(q,w.f\  from  C(0,  T ;  &*)  to  R,  where  /,  is  given 
hi  (4.14)  and  are  the  modal  subapeosa  defined  in  Example  4.2,  is  weH  defined  for 
each  Af  - 1, 2,  •  •  • ,  and,  hi  particular,  is  well  defined  on  solutions  of  (4.12).  Hence 
the  approximating  problems  (ID*)  associated  with  J\  are  well  posed  in  any  event. 
Justification  of  point  evaluation  for  (ID)  with  Jx  depends  heavily  on  the  smoothing 
properties  of  (4.9)  or,  equivalently,  (4.11).  Roughly  speaking,  for  r  >0  die  solution 
values  c(t;«j)  will  be  contained  in  Dom  (A(q))  if  only  t-*f{t)»F{q,  r,  v(t))  is  smooth 
enough  [8,  p.  192].  However,  since  we  wish  to  avoid  additional  smoothness  hypotheses, 
we  choose  a  slightly  more  technical  route  to  the  same  end. 

Recalling  die  arguments  for  Theorem  4.3,  we  have  that  A(q)-ml  is  self-adjoint 
and  maximal  dissipative  where  «  can  be  chosen  independent  of  q  nQ.  It  therefore 
fbHowt  [9,  p.47]  that  -A(«)+wf  mH*{q)  where  «d*(«)  denotes  the  subdifiereatial 
of  the  functions!  d*(q)  is  given  by 

(4.15) 

Here  denotes  dm  square  root,  which  by  standard  results  [21,  p.  281] 

is  known  to  exist.  Under  assumption  (M6**)  for  /,  we  have  that  (H6)  holds  for 
F{q,  U  w)-AV.  s  w)  and  in  particular  (see  (HfiMlH))  foe  mapping  r*F(q,  t,  »(r;  q)) 
is  in  Lg(0,  T\  W)  tor  r-*c(f,q)  the  solution  of  (4.9).  Thai  (see  [9,  p.  72]  end  note  dmt 
a  weak  solution  hi  foe  sense  of  [9]  is  in  fact  dm  unique  mfld  solution  in  our  sense, 
which  is,  moreover,  also  a  strong  sofotion)  h  follows  that  the  map  r-»d*(qKo(r;q)) 
is  in  Li(0,  T;  K)  and  is  foaolutsty  continuous  on  al  subtotarvals  [8,  T\  8>0,  of 
(0,  n  In  |*t  of  (4.13)  fob  kaphas  o(t;q)«Dom(wJ-A(q)),/*  for  all  (>0  and  we 

QOmrmfmm  Wm  mOmm  miMumy  bu  m  m  M,  n  rMnMfil  io  nOtt  mat  Iran  • 
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results  on  interpolating  spaces — see  [27,  p.  9].  Briefly,  by  defining  X  (in  the  notation 
of  [27])  as  either  Hl  with  appropriate  boundary  conditions  (ia  the  event  of  B.C.’s 
using  combinations  from  (A)  and  (B)  or  B.C.(D)  from  (4.13))  or  H1  with  an  appropriate 
energy  norm  (in  the  event  elastic  boundary  conditions  from  (C)  of  (4.13)  are  involved) 
and  Y  *  ft,  one  can  make  the  identification  A  =  (ml  -A(q))l/2  with  Dom  (A)  =  X. 

We  can  therefore  summarize  by  stating  that  conditions  (HQ),  (H4),  (H6**),  (HP1) 
in  the  case  of  (SLBC)  are  sufficient  to  allow  point  evaluation  in  ft. 

If  we  wish  to  relax  the  continuity  hypothesis  on  /  in  Theorem  2.1,  it  is  necessary 
to  consider  in  more  detail  the  fit-to-data  criterion  JN(q)»J(q,  «*  ( • ;  q),  y )  of  (ID*). 
The  following  proposition  will  be  useful  in  our  deliberations;  we  state  and  prove  it 
using  the  notation  of  Theorem  2.1. 

Proposition  4.2.  We  suppose  there  exist  maps/ and /N,  N  =  1,2,  •  • ,  from  the 
compact  set  Q  to  R  satisfying: 

(i)  for  each  N  »  1, 2,  •  •  • ,  the  map  q  -*/N(q)  it  continuous  on  Q-, 

(ii)  for  any  qeQ  and  any  sequence  [Afk]  with  Nk  •*  oo,  there  exists  a  subsequence 
AT*,  such  thatfH*(q)-*J{q); 

(Hi)  for  any  qM  -*4,  there  exists  a  subsequence  (qN‘)  such  that  /tlk(qNk)-*/(4). 

Then  for  each  N  there  exists  qN  e  Q  that  minimizes  /*  over  Q  and,  moreover,  for 
any  convergent  subsequence  {q**}  of  {qs)  with  qMk  -*4,  /  is  a  minimum  over  Oat  4. 

Proof.  Lot  qH  denote  the  minimirer  (whose  existence  is  guaranteed  by  (i)  and 
the  compactness  of  Q)  of/N  so  that/* (4*) 5/* (4)  for  all  4  eQ.  Suppose  qNk  -*#; 
then  by  (iii)  (reindexing  for  convenience  in  notation)  we  have  Jm(qm)-*/(4)  for 
some  subsequence  (4"}  of  {4*}.  We  use  this  to  argue  that  /(4)S/(q)  for  qeQ.  if 
we  assume  teat  there  exist  £  e  Q  and  /  such  that  /m,(qm>)  >/(q)  for  j  2/*  then  by  (fi) 
there  is  yet  anotiier  subsequence  of  mh  denoted  by  mh  such  that/"1"  (9)  -*/(4).  Hence, 
for  sufficiently  large  i  we  have  /m*(qm‘')  >/m“(q),  which  contradicts  the  definition  of 
q"*'1  as  a  minimise!  for 

To  use  Proposition  4.2  with  our  particular  JN  =/N  defined  using  ft,  we  must 
consider  special  cases  for  which  hypotheses  (iHiii)  of  that  proposition  are  readily 
verified.  We  discuss  the  homogeneous  (/—0)  version  of  (4.9)  in  SP  in  this  regard. 

Proposition  4.3.  Consider  (4.9)  wUh  fm0,  initial  conditions  (1C)  and  boundary 
conditions  (SLBQ  from  (4.13).  Assume  dtat  (HQ),  (H4)  and  (HP1)  obtain.  Then  for 
the  solutions  t>w(*;q)  of  the  approximating  equations  ((4.12)  with  F *0)  we  have 
{t>*(q)|q  (O)lta  relatively  compact  subset  of  C(t*,  T;  C(  0, 1;H))  for  each  t*  6(0,  T]. 

Since  the  proof  of  tiiis  result  consists  of  checking  compactness  criteria  for  t  specific 
subset  of  C(t+,  T\  Y),  where  Y  k  a  Banach  spaoe,  we  shall  only  sketch  the  ideas 
involved,  leaving  tire  details  to  the  interested  and  determined  reader. 

First  recall  that  equation  (4.12)  is  in  if’*  -spaa {*1,  •  •  ■ ,  ¥*}  where  1  is 
the  CONS  of  eigenfunctions  of  A(q*)  with  corresponding  eigenvalues  X,  -  Af(q*)  (see 
the  discuaskm  immediately  preceding  Theorem  43).  Solutions  c*  of  (4.12)  have  the 
representation  (see  (2.8M2.12)  and  the  associated  discussions)  c*(f;q)-£*  1  w*(f)y, 
where  wf'(t)  -  w*(0)  exp  {(X<qi+qa)<}-  To  verify  relative  compactness  of  the  desired 
set,  one  can  use  the  Ascoli  theorem  [25,  p.  211]  which  requires  equkontinuhy  of  the 
set  along  with  relative  compactness  of  {vN (f,q)\q a Q,fT -1,2,-  *•,}  in  C(0, 1;  A) 
for  each  t  in  [/*,  T\  Use  of  the  representation  results  along  with  standard  estimates 
in  Fourier  analysis  (Farseval,  Cauchy-Schwsrz,  etc.)  reduces  the  compactness  criteria 
to  the  task  of  verifying 


(W 


1-J 


(c)  -  w*  *  *1,  i  - 1, 2,  *  ■  •}»  flute. 

Here  jf  i.  thetowerboend  for  <j,  (*ee  (HQ)). 

From  tuihdirtf  Fourier  reraht  [IS,  p.  1332]  ee  hive  (A/,  ♦,)-»<)  (or  /€ 
Dom(A(«*))  or  that  k  bOoaried  te  ft  Hence  {♦,}  i>  bounded  In  the 

nri (flit  (Ate*))  norm  end.  u  vnrmwiiiltt  after  tout  aKdMkM  takmt  Into  eccount 
the  (SI-s'C),  In  H\  It  foHow*  that  \c)  bbtafcw  and  that  W}  k  ft*  bounded.  Thu*  if 
(e)  hokb,  We  immedietcty  obtain  (b).  We  have,  therefore,- reduced  (he  compactness 
criteria  to  a  requirement  oo  the  rate  at  which  X,  -*-eo.  Hot  aiymptodc eethnatee for 
toe  dfenvaiues  or  sctBiiHJoiivuR  pfoown  rcooHy  svatisDw.  roc  csnnpfCy  tor 
(SLBO  from  combination*  of  (A),  (B),  (Q  of  (4.13)  we  have  h — -/*  [14,  p.  1531 
which  yield*  (a).  For  the  periodic  boundary  peufltfooi  (D)  a  slight  modification  of  the 
aifltbfliftrftt  [12,  §  8.3]yiekJ  the  «ame  rate  edflmat**.  In  aome  higher -dimensional 
problems,  simitar  rates  for  the  eigenvalues  are  aho  evnflahlc  [13,  *  VL4.1]. 
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4.3  we  ate  npw  able  to  give  well- 
0  fa  (4.9)- (IC)- 


4-*  I  fC(44)«ol(ew(4*i;f),*‘-,ew(lW*t;f)>-,r(4# 
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•  ■■  ■ 

fi*  4**  with  4K+4*  Q 

HtMlstharwfca  *ifltTarilg  shown  hi  (Q,  T),  we 
‘Mhem^ewee  {»**»(•  ;4W‘)) 


wOethy  in<C«%2>,  C(Q.  !;&)).' 


)-4<tt4)fc^«dl*a#inron[#*,  T\ 

of  the 
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Mflifft**  biebi!*;!?  dire  i.-rit* 
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the  operator  A(f)  by 
>4**.Asnah«(HQ)aad<H4>. 
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it  is  easy  to  verify  existence  of  an  <u  independent  of  q  e  Q  such  that  A(q)-a>I  is 
dissipative  and  symmetric  in  that  Dom  ( A[q ))  is  dense  and  R(A(q)-kI)  =  &  for 
some  A,  chosen  sufficiently  large  and  independent  of  q  €  Q.  In  particular,  A{q)  is 
therefore  self-adjoint  [28,  p.  97].  The  approximating  subspaces  are  defined  by 

-  span  {4>i,  •  ■  • ,  where  4>/(x)  *  (v2//w)  sin  jvx  as  in  Example  4.1.  We  recall 

that  is  a  CONS  in  &.  As  before  we  let  denote  the  canonical 

orthogonal  projections  and  AN(q)**PNA(q)PN,  FN(q,t,v)**PNF(q,t,v)  where 
F:  Q  x  [0,  T]  x  %•*  £  is  given  by  F{q,  t,  c)=  /(<?,  t,  • ,  o).  Formulation  in  X  will  require 
additional  assumptions  on  /,  to  be  detailed  in  (H6***)  below.  We  then  have: 

Theorem  4.5.  Suppose  that  (HQ)  holds  and  let  qs,ijeQ  be  such  that  qN-*q 
as  IV -*oo.  Then  A (4)  and  AN(qN)  generate  Co-semigroups  T(t;  q)  and  TN(t; qN)  on 
&  and  \TN(t;qN)z  —  T(t;4)z|jr-»0  for  each  ze&,  with  the  limit  uniform  in  t  on 
compact  subsets  of  [0,  oo).  Furthermore ,  if  Q  is  compact,  then  there  exists  a  constant  o> 
independent  of  N  and  q,  such  that  |r(f;q)|Se“'  and  \TN(f,q)\SeM  for  t>0,  qeQ, 
and  AT  *  1, 2,  •  •  • . 

The  proof  of  this  theorem  is  quite  similar  to  that  of  Theorem  4.3  and  will  therefore 
not  be  given  here. 

We  shall  simply  list  conditions  (compare  with  (H6**))  on  /  that  will  insure  that 
F  satisfies  (H6). 

(H6***)  The  nonlinear  function /:  O* * [0,  T]x[0,  l]x R-* R  satisfies: 

(0)  For  each  v  €  the  map  x  -*f(q,  t,  x,  u(x))  is  in  ffc 

(i)  For  each  (q4,t>)€Q4xR,  dm  maps  (r,x)-*/(q*,t,x,i>),  (t,x)-* 

fn(q*.  t,  x,  o),  and  (/.  x)-*f,(q*>  <.  *,  »)  are  measurable. 

(ti)  There  exists  a  function  k\  in  L2(0,  T )  such  that  \fdq*<  A  x,  v)ls£t(t) 
for  all  (q4,  x,  o)e  Q4x  [0,  l]x/?;  for  each  At  >0  there  exists  in  L2( 0, 1) 
sudi  that  \fx(q*,  t,  x,  vi)-fx(qt,  t,  x,  o2>!  *£i(t)\vx  -t>2|  and  IfAq*,  t,  x,  ui)- 
Mq*.  t,  x,  t>2)|  3  k*(t)\vi- u2|  for  all  (q4,x)e<?4x[0, 1]  and  ci,  o2  with 
|t>i|SAf,  |c2|SAf. 

(iii)  There  exists  k2  in  L2(0,  T)  such  that  |/(q«,  r,x,0)|S*2(f)  and 
l/«(q«>  ^ x,  o)|  S  £2(t){l + |u|}  for  all  (q4,  x,  i>)eQ4x[0,  l]xjf. 

(iv)  For  eadi  (t, x,  v)  in  [0,  r]x[0,  l]x/f,  the  maps  q4-f{q4,  t, x,  v),  q*-* 
fxlq*,  t,  x,  o)  and  q*-*f,(q*,  t,  x,  v)  are  continuous  on  Q4. 

The  ftt-to-data  criterion  (4.14)  together  with  the  state  equation  in  #  are  such 
that  the  map  (q,  v)-*/i(q,  v,  f)  from  Q  x  C(0,  T ;  &)-*R  is  continuous  and,  therefore, 
Theorems  3.1  and  2.1  are  readily  applicable.  We  leave  a  precise  statement  of  the 
theorem  for  (4. 16MICHDBC)  with  (H6***),  analogous  to  Thoorem  4.4,  to  the  reader. 

in  concluding  this  discussion,  we  remark  that  numerical  implementation  of  a 
scheme  formulated  as  above  in  4  (the  projections  in  defining  the  approximations  in 
(2.8H2.12)  are  now  in  the  o  inner  product)  is,  of  course,  somewhat  more  tedious 
from  a  technical  viewpoint  than  that  for  schemes  such  as  those  in  Example  4.3  where 
the  state  spaee  is  H°. 

Example  4.4.  Diffur'on-convection  equations.  For  the  final  example  of  this  sec¬ 
tion,  we  return  to  the  setting  of  Example  4.2  and  indicate  how,  in  (4.9),  one  might 
indude  ocmvoetlon  (or  advection)  terms  that  are  independent  of  the  Sturm-Liouville 
operator  (pa,)*  To  illustrate  toe  ideas,  we,  for  ease  in  exposition  only,  take  a  simple 
Unear  example  (nonlinearities  of  toe  type  discussed  previously  present  no  essential 
dlfficuHlss)  involving  only  diffusion  and  convection  terms.  Consider  then 


(4.17) 
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with  initial  conditions  (IC)  and  the  Dirichlet  boundary  conditions 
(DBC)  t>(f,0)  =  t>(f,l)  =  0. 

As  state  space  we  choose  H°  with  its  usual  inner  product.  We  define  the  operator 
A(q)  with  Dom  (A(q))  =  H2  Cl  Hi  by  A(q)ifr  =qtif'”+q2</''-  We  then  have  A(q)  is 
dissipative  since  (again  we  assume  (HQ)  holds)  for  z  e  Dom  (A  (4)), 

(A(q)z,  z)  =  qi(zxx,  z)+qi(zz,  z) 

=  -<7i|z*|2+<f2  [  (z2/l)xdx 
Jo 

*-^i|r,|2S0. 


Thus  if  we  assume  (H4),  we  find  that  A(q)  is  dissipative  uniformly  in  q  e  <?. 

In  considering  modal  approximations,  the  question  erf  existence  of  a  complete 
set  of  eigenfunctions  for  the  operator  A(q)  arises  naturally.  Standard  spectral  results 
for  nonself-adjoint  operators  allow  one  to  answer  this  question  in  the  affirmative. 
First,  A(q)  is  a  relatively  bounded  perturbation  of  a  discrete  spectral  operator  and  is 
itself  a  discrete  spectral  operator  (see  [IS,  Thm.  XDC4.16,  p.  2347] — in  this  case  the 
boundary  conditions  (DBC)  are  easily  seen  to  satisfy  the  necessary  regularity 
hypotheses — see  [IS,  p.  2341-2344]).  It  follows  that  <r(A(q))  consists  of  point  spectrum 
and  that  the  eigenprojections  {£*,}  (see  [IS,  p.  2292])  of  the  resolution  of  identity  for 
A(q)  satisfy  if.j  E^z  -*z  for  z  eH°  (see  [15,  Cor.  XVffl.2.33,  p.  2257],  along  with 
the  properties  of  the  projection  operators— e.g.,  [IS,  Lem.  XVIU.2.31,  p.  2255]). 
One  can  easily  argue  for  our  example  that  the  generalized  eigenmanifolds  are  one- 
dimensional  so  that  the  eigenfunctions  $;(q)  =  exp  (-q&flqi)  sin  jirx  corresponding 
to  the  eigenvalues  A t(q)  =  -j2ir2qt -ql/2qx  form  a  complete  (but  not  orthogonal)  set 
in  H°. 

We  thus  also  have  that  A  ep(A(q))  if  A  >0  so  that  (A(q)-l)  Dom  (A(q))  =  H° 
for  A  >  0  and  hence  A(q)  is  maximal  dissipative  [22,  p.  87],  [30,  p.  17].  The  operator 
A(q)  generates  a  Co-semigroup  T{t; q)  satisfying  |7'(<;q)|Se“'  for  qeQ. 

For  a  modal  approximation  scheme,  it  might  be  tempting  at  first  thought  to  use 
the  finite-dimensional  subspaces  £N(q) «  span  {♦1(9),  •  •  • ,  ^n(9)}>  but  of  course  this 
would  prove  rather  difficult  computationally  in  identification  problems.  Here  we  choose 
to  use  the  basis  elements  ♦/(x)»>/2sin/irx  since  we  know  forms  a  CONS  in 
H*  and  Q/eDom  (A(f».  We  thus  define  HN  -  span  {4>i,  •  •  • ,  <t>v}  and  remind  die 
reader  Oat  “modal”  is  something  of  a  misnomer  for  this  scheme  (actually,  we  took 
a  shnBar  approach  in  Example  4.1  hi  choosing  basis  elements  corresponding  to 
*•-(1,0,  •  --.O)  fixed). 

As  usual,  we  define  AM(q)  -  ^AWP"  where  P*  are  the  orthogonal  projectors 
PNz  -Ejij  (x,  ♦/)♦/  onto  H"  which  converge  strongly  to  the  identity  on  H°. 

To  develop  approahnatfoa  results  shatter  to  those  given  in  Theorems  4.3  and 
4.4,  (ha  eamattal  efloct  rsasaWag  in  oer  enampla  is  to  verify  the  stability  and 
cooaMswcy  hypotheses  (R),  (if),  (hr)  of  Proposition  4.1  with  j rf- A(«f)  and  siN  - 
A* (4*)  where  9*  in  Q.  StaMky  (Le.,  (if))  is  Immediate  while  consistency  is  slightly 
more  dottcasa.  A  animal  Asin  Mm  one  we  have  used  in  nrevioBs  examoles)  for  the 
set  9  is  U  7.|  MN  Mam  than  0v)  is  trirtf  to  verify.  However,  it  is  not  apparent  to 
as  that  tMs  choke  of  9  satkdss  (id)  of  PropastHiin  4.1.  We  choose  inmead  9  - 
UU.i^d),  where  £*<#)  h  daAnod  above  in  terms  of  dm  true  modes  +,(4)  for 
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A(q).  Since  (Al -A(4))4'/(4)  =  (A -A,)^,(4)  and  the  set  {^(4)}  is  complete,  (iii) 
is  easily  established  and  it  only  remains  to  show  AN(qN)z  -*A(4)z  for  z  e  2. 

We  first  note  that  from 

I AN(qN)z -A(q)z\H\PN(A(qN)PNz -A(q)z)\  +  \PNA(4)z -A(4)z\ 

and  the  strong  convergence  of  PN  to  /,  it  is  sufficient  to  argue  A{qN)PNz  -*A(q)z 
for  ze2.lt  suffices  to  argue  this  latter  convergence  for  z  =  fa  »'i'*(4)  fixed.  For  this 
choice  of  z  we  find 


AN(qN)PN*k=A(qN)  £  <^,4*^=  £  {^,)A(qN)<bi 

i- i  /-i 

=  z 

/-i 

=  I  +iirqzXi} 

l-l 

=  qi  Z  (4k,-j2ir2$i)Qi+q2  Z  (<l>k,i7r<t>,)xi 

I- t  /-i 

=41*  Z  (dfe  4»/')4>/  +4*  z  (4>k,  -x))Xb 

where  *y(x)  =  V2  cos/jtx  and  we  have  used  the  facts  that  x/=-/V4>>  and  4>*  = 

-/V*,. 

Integration  by  parts  twice  (using  the  fact  that  &  and  4>,  are  in  //o)  yields 

4>'>  =  (</>*, 

while  a  single  integration  by  parts  establishes  (again  use  e  Ho) 

Wk,-x’l)  =  ('l'k,Xi)- 

We  thus  have 


AN(q")Ps*k-q?  1  <#S,  Z 

/-i  /-i 


where  both  {♦/}  and  {*/}  constitute  CONS  in  H°.  Since  q*  -*4u  q? -* 4z  we  thus 
obtain  AN(qt,)PNif>k-*4nl>k  +4j*»  -  A  (4)+*,  as  was  desired. 

The  theorems  for  these  approximation  ideas  for  the  (ID)  and  (IDN)  problems 
with  (4.17HICHDBQ  are  so  similar  in  statement  to  Theorems  4.3  and  4.4  that  we 
shall  not  prolong  our  discussion  by  giving  a  precise  statement  here. 

With  regard  to  implementation  of  this  scheme,  we  point  out  that  A(g)  does  not 
leave  the  subspacesff11'  invariant  and  hence  the  matrix  representation  of  A*  -  PNAPr> 
(see  (2.8H2.10))  is  not  a  simple  diagonal  matrix.  However,  for  equations  such  as 
(4.17),  it  is  rather  easily  seen  that  (2.10)  is  given  by 


£*"(«>]*- 


for/»/, 

for/  */ and/  +/  even, 
for/**/ and/+/  odd. 


While  this  is  not  a  staple  matrix,  it  does  aQow  a  rather  straightforward  implementation 
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5.  A  boundary  control  probka.  The  theory  developed  in  §§  2  and  3  can  also 
be  applied  to  optimal  control  problems  governed  by  partial  differential  equations.  We 
shall  demonstrate  this  by  means  of  a  specific  example.  Consider  as  a  special  case  of 
(4.1HICiMBCi)  the  problem 

(5.1)  va  =  v„ 

for  t  >  0,  x  e  [0, 1]  with  initial  and  boundary  conditions 

(ICj)  J(0,  x)~4>{x),  v,(0,x)  =  i(i(x), 

(BC3)  v(t,  0)  =  v(t,  1)  =  j2(f), 

where  the  boundary  control  functions  s(  are  chosen  in  5^={s|s6/f2(0,  T;R), 
s(0)  —  j'(0)  =  0}  and  (<t>,  4>)  e  Ho  x  H°.  The  transformation  v(t,  x)  = 
v(t,  x)-(l -x)si(r)-xs2(f)  applied  to  (5.1HIC3)-(BC3)  leads  to 


(5.2)  »«  =  »*, -(l-x)(*i)«-x(*2)«, 

(IC4)  u(0,x)  =  ^(x),  u,(0,x)  =  dr(x), 

(BC4)  o(t,0)*»(r,  1)  =  0. 

We  let  w  =v,  and  reformulate  (5.2)-(IC«HBCi)  as  in  Example  4.1  in  the  Hilbert 
space  X=Hh*H°  with  the  usual  inner  product.  This  leads  to  a  special  case  of  (4.3) 
given  by 


where 


£),  y-y(x)-(_°  <r(r)«co!  ((,,)«.(,,)„).  (d,if>)eX. 

The  finite-dimensional  subspaces  Xs  »  Xs (q*),  q*  =  (1, 0,  •  •  • ,  0),  are  chosen  as  in 
Example  4.1  and  again  we  take  A  w  =»  PNAPN,  where  PN :  X  -*  XN  denote  the  canonical 
orthogonal  projections.  For  the  convenience  of  the  reader  we  repeat  the  family  of 
approximating  equations  given  by 


(5.4) 


In  the  light  of  Theorems  3.2  and  4.1  (with  qM  mq*  for  aOM),  the  solutions  (c*(f;0-), 
w"(t;  <r))  and  (c(f;  or),  w«;  <r))  of  (5.4)  and  (5.3),  respectively,  satisfy  tim*  (cH(t;  <r), 
w"(t;  <r))  -  (»(»;  <r),  w(f;  <r))  in  X  unfformly  in  t  «R),  T],  for  any  T>0  and  uniformly 
in  <r,  as  o-  varies  in  bounded  subsets  X  of  L2(0,  T\  R1).  We  shall  also  need  the  foilowii^ 
technical  result. 

Lemma  5.1.  The  operator  XiLdO,  T,Ra)-C( 0,  T;X)  defined  via  (ArKO* 
£  T(t-r)ytr(.r)  dr  is  compact. 

Defining  the  maps  ■  £  T*(t  -  r)FNyo(r)  dr  and  using  the  convergence 

of  the  semigroups  T*{t)  to  r(r),  generated  by  A "  and  A,  respectively,  it  is  easily 
seen  that  hi  the  operator  norm  topology.  The  proof  is  completed  once  one 

argnse  that  the  maps  X"  themselves  are  compact 

IN  WWfJ  IWBUU  psvviw!  IBf  Wflilwll  wOH  UNU  Oil  M  WppaBQ  ID  ■  tiiigij 

of  optimal  control  problems,  ooc  of  which  vil  be  otttHned  below.  Bor  a  more  complete 
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discussion  concerning  approximation  of  optimal  control  problems  for  infinite¬ 
dimensional  systems  by  sequences  of  optimization  problems  for  finite-dimensional 
systems,  we  refer  to  [7]  and  the  references  given  there. 

We  let  £«)  be  any  fixed  dosed  convex  subset  of  L2( 0,  T\ R1)  (possibly  L2  itself) 
and  choose  nonnegative  continuous  functions  g0:9t-*R,gi:  C(0,  T ;  X)-*R  and 
ft: £2(0,  T ;  /?*)-*/?.  The  functions  g(  define  the  cost  functional  /:  X*) -*R  by 

(5.5)  S(<r) » gMp(T ;  a),  w(T ;  er)))+gi((e(  • ;  tr),  tv( • ;  o-)))  +  g2(o). 

The  optimal  boundary  value  control  problem  assodated  with  (S.lHICaHBCj),  (5.5) 
is  then  taken  to  be: 

(iP)  minimize/ over 

Suppose  that  a  solution  <r = col  (di,  d2)  6  of  (9)  is  found;  this  will  uniquely 
determine  boundary  controls  fi  and  f2  in  S'*.  The  approximate  optimization  problems 
are  defined  by 

(&N)  minimize  JN  over 

where 

(5.6)  /*(*)- go((»N(T;o-),  w"{T\  <r)))+g,((eN(- ;  w),  *vN(- ;  «•)))+ *2(«r). 

Notice  that  (5.6)  is  an  optimization  problem  associated  with  an  ordinary  differential 
equation.  We  shall  need  two  standard  assumptions  on  die  functions  g,i 

(Ol)  The  continuous  functions  g<  are  convex, 

(GT)  gi(o)  -*oo  as  t<rj-*oo. 

As  a  consequence  of  (Ol),  the  maps  cr-*/(er)  and  <r-*SN(<r)  are  convex,  which 
together  with  (02)  implies  the  existence  of  solutions  of  (P)  and  (PN);  these  solutions 
are  in  addition  unique  if  one  of  the  g<  is  strictly  convex.  Let  <rN  denote  a  solution  of 
(^N).  Then  by  (02)  it  follows  that  {o-*}  must  be  a  bounded  subset  of  L2(0,  T:  R2). 
Indeed,  the  assumption  jo^l-vao  for  some  subsequence  {Nk}  contradicts  die 
inequalities  g2(«rw*)a/Af,‘(orw*)JS/Wk(cr)-p/(o,)<«>  for  all  <r*Lj(0,  T;R2).  The 
convergence  of  JNh(<r)-*J(<r)  is  a  consequence  of  (»w(r;cr),  ww(f;<r))-* 
(»(/;<r),  w(t;a))  uniformly  in  ie[0,  T]  and  (Ol).  Since  X*  is  convex  and  closed  it 
is  weakly  closed  so  that  there  exists  a  weakly  convergent  subsequence  {?”*}  of  {<rN) 
with  <rN*  converging  weakly  to  some  d  e  W  By  Theorem  3.2,  Lemma  5.1  and  the 
estimates 

|(oN<;  <r*‘),  wN‘(f;  <r'f‘))-(p(t ;  d),  *v(»;  d))| 

S|(e N»(r;<rw*),  wN*(t\<rNk))~(v{.f, <rN“),  w(l; crw*))| 

+|(t>(r;<rw»),  w(t\<rSk))-(v(t\&),  w(f;d))|, 

it  follows  that 

(cN‘(»; <r"'),  ww*(r;  tr"*))-* (»(r;  d),  w(f; d)) 

in  5T  uniformly  in  r  a  [0,  T).to>ce  convexity  and  continuity  together  imply  weak  lower 
semicontinuity,  we  obtain  the  fbOowing  string  of  inequalities: 

/(d)*hmtaf  ige((»w*(r;<rw*),  ww*(T;<rw*))) 

♦#!«»*(•  i  w"*),  *>*(•;  w*)»+ffa(w"k)} 

-  Um  faff  /Nw1 N»)  a  Ua  sup  /"*<*)  ■  /(w) 
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for  every  <r  e  San.  This  implies  that  &  is  a  solution  of  (fF).  Further  standard  arguments 
can  be  we d  to  show  that  strict  convexity  of  J  implies  dim  aN  itself  converges  weakly 
to  the  unique  solution  <r  of  (9)  an d  that  <rN  converges  strongly  in  Lj(0,  T\  l!2)  to  a 
if  /  is  strongly  convex  (see  [7]).  We  finally  summarize  some  of  the  above  discussion  in 

Theorem  5.1.  Suppose  that  (Ql)  and  (G2)  held.  If  {<?*}  denotes  a  sequence  of 
solutions  of  (9N),  then  there  exists  a  subsequence  {<rN‘}  converging  weakly  to  a  solution 
d  of  (9).  Furthermore,  JNk(<rSk)-*S{&)  and  (vr,k(tnrflk),wNk(t;<rNk))-* 
(t >(t;o-),w(t;&))  uniformly  in  fe[0,  T].  Moreover,  &  determines  uniquely  boundary 
controls  S\,  i3  in  9. 

6.  Numerical  rramplts.  In  this  section  we  briefly  summarize  our  numerical 
findings  when  applying  the  modal  approximation  algorithms  to  some  of  the 
identification  problems  that  were  outlined  in  8  4.  The  aim  here  is  to  demonstrate  the 
feasibility  of  die  method  for  both  hyperbolic  and  parabolic  systems.  As  it  turns  out, 
modal  approximations  appear  to  be  very  well  suited  for  hyperbolic  systems,  while  for 
certain  identification  problems  for  parabolic  systems  we  encountered  some  essential 
difficulties  which  one  should  take  into  consideration  before  attempting  any  practical 
use  of  the  method  for  this  type  of  equation.  This  will  be  explained  further  below.  In 
developing  our  software  packages,  no  great  attention  was  given  to  maximizing 
efficiency  in  implementing  the  algorithms,  or  to  minimizing  computer  time.  The 
ordinary  differential  equations  (see  (2.12))  that  arise  were  integrated  by  a  simple 
fourth-order  Runge-Kutta  method  (with  step  size  varying  from  one  example  to  the 
next  from  .0125  to  .05),  and  the  coefficients  of  the  nonlinearity  and  the  initial  data 
(see  (2.9)  and  (2.11))  were  computed  by  employing  Simpson’s  rule.  The  minimization 
problem  arising  in  the  identification  problem  for  the  approximating  ordinary  differen¬ 
tial  equations  was  numerically  solved  by  using  an  IMSL  package  (ZXSSQ)  employing 
the  Levenberg-Marquardt  algorithm.  The  “exact”  solutions,  which  were  used  for  the 
“data"  £  in  the  fit-to-data  criterion  /,  were  generated  by  a  Crank-Nicolaon  algorithm 
whenever  solutions  in  closed  form  were  not  available.  These  solutions  were  generated 
with  fixed  known  values  of  the  parameters  in  the  equations;  these  values  will  be 
referred  to  in  the  sequel  as  the  "true"  parameter  values. 

In  the  examples  below,  a  fit-to-data  criterion  of  the  type  (2.4)  with  C(thq)~I 
was  used  throughout.  Further,  we  usually  (except  in  Example  6.5)  let  T  ■  2  and  chose 
t,  and  x;  equally  spaced  in  [0, 2)  and  [0, 1],  respectively,  so  that  =  0.2  and 

|x,-x,_,|-0.25. 

Example  6.1.  Here  we  return  to  Example  4.1  and  consider  the  linear  one- 
dimensional  hyperbolic  equation,  which  we  repeat  for  convenience: 

n**9i»*»+q3U,+qjt;  for  f  >0, 

t>(0,x)»4«x(l-x)  forOSxSl, 

c,(0,x)-qsd(x)  forOSxSl, 

c(t,0)-e(t,  1)*0,  forr>0, 

where  tf(x)-2x  for  xc[0,  .5]  and  <f(x)-2(l-x)  tor  x  c  [.5, 1J.  Below,  we  present 
numerical  results  which  aw  typical  of  those  obtained  hi  making  numerous  runs  with 
this  example.  The  startup  values  qm  for  use  hi  the  Levenberg-Marquardt  algorithm 
an  recorded  in  the  bottom  row  of  the  tables,  whereas  the  next-to-last  row  contains 
the  true  parameter  values.  Urn  tables  contain  only  those  parameters  on  which  a  search 
was  performed,  whewas  the  wmahriag  parameters  wew  assumed  known  and  therefore 
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In  a  first  run  (Table  1)  we  assumed  that  q3  =  1  and  q3  =  0  were  known  and  we 
were  searching  for  q  =  (41,  q2,  q4)  with  the  true  parameter  values  chosen  to  be  4  ~ 
(1.414,  -1, 2)  am)  with  startup  values  qN'°  =  (1, 0, 1). 


Table  1 


iV 

4x 

4 r 

4 

1.4103 

-0.9961 

1.9978 

8 

1.4126 

-1.0021 

2X1031 

16 

1.4129 

-0.9968 

2.000S 

32 

1.4129 

-0.9992 

2.0000 

true  value 

1.414 

-1 

2 

Is* 

1 

0 

1 

A  feature  of  interest  for  these  models  used  with  die  Levenberg-Marquardt 
algorithm  is  the  range  of  convergence  for  the  parameter  q.  For  this  specific  example, 
we  carried  out  computations  keeping  two  of  the  parameters  41,  42,  4«  in  addition  to 
43, 4s  fixed  while  identifying  one  of  41, 02, 44-  It  was  observed  diet  for  <jf°,  4?’0, 4  J'  0 
taken  in  the  ranges  1  Sq%  S3,  -SS42  SO,  .5  SqJ^  *3,  respectively,  rapid  conver¬ 
gence  was  still  obtained.  (The  actual  range  of  convergence  may  be  much  larger;  these 
were  merely  the  ranges  of  values  sse  tested.) 

In  a  second  run  (Table  2),  41  - 1.414  mid  4*- 1  were  assumed  to  be  known  and 
the  search  was  performed  on  4  •  (42. 43, 44)  with  true  value  4  “  (-3, 4, 2)  and  startup 
values  q"fi  *  (0, 0, 1). 


Table  2 


N 

4? 

4? 

4? 

4 

-4.9*30 

4.0342 

1.99*7 

S 

-4.9003 

4.0972 

2.002S 

trasvslas 

-s 

4 

2 

9*** 

0 

0 

1 

Example  6.2.  This  is  the  aoaMnaar  example  (^Sa  a  pdd  case  of  Busmgls  4.1): 

e(0,x)-g4*(l-x)  •arOSxSl, 

c,(0,x)-0  farOSxSl, 

ett0)-»(fcl)-0  lerr  >0. 

We  chase  the  true  model  pemesassss  4  -  (1.414, 2, 1),  whereas  dm  startup  values 
were  taken  to  be  4**  -  (1, 1, 0).  For  dss  usmirlssl  tahsiens  we  rslsr  to  Tshia  3. 
Bxampk  6.3.  TMsis  sne*ar  nauMnaar  efuatfanef  thefem 

as -flew 4-4,0 «►#•»*  tsrf>0» 

»f0,x)-44«(l-*)  hfOligl, 

***)«•  lurOSxSl, 

*&•)-*&  l)-0  tsrf  >0. 


i 
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Table  3 


N 

4 1 

4? 

4? 

4 

1.4141 

1.9990 

0.9733 

8 

1.4148 

2.0013 

0.9790 

16 

M1S2 

2.0009 

0.9788 

true  valve 

1.414 

2 

1 

qK* 

1 

1 

0 

Although  this  aonhaeshty  does  not  satisfy  (H6*)  of  Example  4.1,  we  report  in 
Table  4  on  calculations  carried  out  in  the  subspaces  XN  of  X.  It  is  clear  that  die 
algorithm  is  converging  in  this  case;  indeed,  one  can  relax  the  assumptions  (H6*)  so 
as  to  prove  convergence  for  such  nonlinearities;  see  the  discussion  involving  (A6)(i), 
(ii)  in  S  3. 

Table  4 


N 

4t 

4j 

4? 

4* 

4 

1.3833 

0.6774 

1.9999 

1.2368 

8 

1.4107 

0.987S 

2.0001 

0.8973 

16 

1.4138 

0.9983 

2.0001 

1.0016 

tnw  value 

1.414 

1 

2 

1 

9 

1 

0 

1 

0 

We  turn  now  to  some  special  cases  of  die  parabolic  problem  (4.9MICHBC).  As 
pointed  out  earlier,  parabolic  equations  can  be  more  formidable  than  hyperbolic  ones 
to  handle  via  modal  approximations.  The  difficulties  are  more  than  just  a  simple  lack 
of  idendfiability  (however  this  concept  is  defined),  which,  of  course,  can  lead  to 
substantial  numerical  embarrassment.  Indeed,  parabolic  equations  can  lead  to  stiff 
systems  of  approximating  ordinary  differential  equations.  The  reader  can  quickly 
convince  himself  of  this  fact  by  taking  Dirichlet  boundary  conditions  and  putting 
and  fmqtmO-  In  our  computational  pursuits  we  did  not  make  an  effort  to 
use  specific  numerical  methods  for  the  stiff  systems  that  can  arise,  but  we  simply 
decreased  die  step  rise  in  die  Ruage-Kutta  algorithm  to  effect  numerical  stability.  A 
perhaps  asore  reasonable  approach  to  avoiding  these  difficulties  due  to  modal  approxi¬ 
mations  is  to  take  a  completely  different  approximation  scheme,  say  for  example 
spliae-based  smthods.  We  have  pursued  tide  idea  tuceessfuffy  for  parabolic  systems 
and  the  details  of  those  investigations  wifi  be  reported  elsewhere. 

The  flt-to-data  criterion  is  chosen  to  be  (4,14)  with  C(rhq)«/  in  all  the  scalar 
■ |«t—  below.  In  the  two-dimensional  system  of  Example  6.7,  we  used  the  obvious 
analogue  of  (4.14)  for  a  coupled  system  of  equations. 

Exampk  6.4.  We  consider  the  linear  equation 

<’fa>4iBtt+42e  forr>0, 

u(0,x)-^C*)  forOSxSt, 

p(A0)-e(/,  l)-0  forr>0, 

where  f  is  the  "hat” -function  defined  in  Example  6.1.  The  modal  approximation 
scheme  faffed  to  identify  ft  and  43  limuhensnusly,  although  it  did  identify  each  of 
them  todlridustty  so  long  m  the  other  one  was  fixed.  Tbh  is  by  no  means  surprising; 
dm  exact  solution  of  the  above  problem  hat  the  expffdt  representation  e(t,x)« 
X£i  c,  (r)  tin /ex,  where  uf(r)»u,(0)eBp((f|-ft(/»)*)<)end  c,(0),/- 1,2,  •  •  •  .are  the 
Rentier  coefficient!  of  the  tine  aariee  for  £  At  thne  t  -  0,  the  values  of  qu  qt  have 
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no  influence  and  at  the  following  times  *  -  .2,  .4,  •  -  the  decay  of  the  exponential 
term  in  addition  to  the  decreasing  magnitude  of  the  coefficient  e/(0)  erase  successive 
terms  to  contribute  less  to  the  criterion  /.  Moreover,  in  this  example,  o2/(0)-0  for 
/  *  1,  2,  •  •  • ,  so  that  the  criterion  uses  essentially  only  one  mode  to  (it  the  model  to 
the  data.  The  results  from  the  search  on  both  parameters  simultaneously  are  presented 
in  Table  5. 

Table  S 


N 

4" 

4 7 

4 

0.0236 

0.2313 

8 

0.033S 

0.3289 

16 

0.0336 

0.3296 

true  value 

0.1 

0.986 

0.23  . 

0.2S 

Keeping  o2  -  .2  fixed  and  searching  for  qu  when  die  true  parameter  value  is 
4i-.l  and  4i  -  0.2S,  we  find  -  0.09999.  Similarly,  when  41  ->  .1  is  kept  fixed  and 
42  is  to  be  identified,  with  “  -986  and  q"'°  -  .23,  the  algorithm  yields  986004. 
Example  6.S.  We  next  consider  die  nonlinear  parabolic  equation 

v, -qiv„  -44e*  for  r>0, 

v(0,x)-q3if  forOSxSl, 

o(t,0)“»(t,  1)“>0  fort>0. 

It  is  well  known  that  for  4«k0  the  above  system  has  a  global  solution  and  we  are 
therefore  again  in  a  situation  where  hypotheses  (A6)(i),  (ii)  of  1 3  must  be  used  in 
any  theoretical  considerations  of  convergence.  Our  findings  for  dds  example  are  given 
in  Table  6.  Here  we  choose  T  -  1,  while  keeping  die  increments  between  die  “data” 
points  the  same  as  before  (A t « .2  rad  Ax  -  .25). 


Table  6 


N 

4i 

4? 

47 

2 

.3030 

4.8271 

.8339 

4 

4976 

5. 3001 

1.2370 

a 

.4983 

3.1774 

1.1482 

16 

.3021 

3.0843 

1.0443 

true  rah— 

.3 

5 

1 

25 

1 

0 

Example 


below,  •  l 


6.6.  We  consider 

»."4i»*«+2q«(l+»rt  fort>0, 
c(0,  x)«*4j#  forOffixffil, 

for  r>0. 

we  solved  the  approximating  identification  problem  both 
added,  then  the  Graak-Nteohao  data  which 
criterion  were  perturbed  by  Gaussian  noise  with  aero 
ee  #*-.01.  It  is  accurate  to  report  that  in  these  two  examples  die 
in  a  stable  maimer  under  the  influence  of  noise.  In  Tabiee  7  and  8 

was  kept  find  .at  the  true  parameter 


f 


-  •  .  -  m"*"'  '  ’ 

— fttea  a 


r _ 

te.-ssaPfi 


mam.:** 
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Tablb  7 


N 

4? 

4? 

4? 

no  note,  itudi  oo  q},  q4 

4 

S.227S 

1.9254 

8 

5.1374 

1.9741 

16 

S.0668 

1.9845 

note  cr’-.Ol;  search  on  4],  q. 

4 

S.2362 

1.9335 

S 

5.1439 

1.9813 

16 

5.0749 

1.9917 

no  note;  search  on  41,  4j,  4< 

4 

0.2472 

5.2846 

2.5221 

8 

0.2301 

5.1706 

2.3584 

16 

0.21  SO 

5.0823 

2.1746 

note  <r2-.01;  search  on 41, 41, 44 

4 

0.2443 

5.2903 

2.4941 

0.2272 

5.1760 

2.3301 

16 

0.2120 

5.0873 

2.1442 

trae  vitae 

0.2 

5 

2 

KO 

4 

0.1 

1 

0 

Example  6.7.  As  a  final  example  we  consider  the  coupled  parabolic  system 

v,  +2(1  +q4w+v)'\ 

for  f  >  0, 

e(0,  *)-#(*) 

forOSx  SI, 

w(0,x)-^(x) 

forOSx  SI, 

0(1,0)  -0(41)- 

w(l,0)»w(t,  l)-0 

forr>0, 

lor  which  the  numerical  results  are  given  in  Table  8. 

Table  8 

N 

4? 

4* 

no  none;  tearch  on  41, 44 

4 

.2011 

1.9933 

8 

.1087 

2.0226 

note?*  *.01;  search  on  41,4* 

4 

.1982 

2.1103 

8 

.1960 

2.1246 

no  note;  search  on  43, 44 

4 

.0500 

2.0514 

8 

.0498 

1.9551 

note  r*«,01;  search  00  4*  «« 

4 

.0522 

2.0349 

8 

.0520 

1.9385 

no  note;  search  oa  43, 43,  q4 

4 

.2011 

.0500 

1.9931 

• 

.1988 

.0499 

2.0187 

ante  »*“.01:  search  00  43,  ft.  ft 

4 

.1973 

.0522 

2.1776 

• 

.1949 

.6521 

2.2066 

tmi  viloft 

.2 

.05 

2 

^MS 

.1 

.1 

0 
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7.  Cendodlag  wta  The  contributions  of  the  discussions  in  this  paper  are 
twofold.  First,  we  have  developed  a  general  approximation  framework  in  the  context 
of  semigroups  that  allows  treatment  erf  identification  and  control  problems  for  a  wide 
class  of  distributed  parameter  systems.  Our  second  contribution  is  a  development, 
using  this  framework,  of  “modal"  approximation  schemes  in  the  spirit  of  those  often 
proposed  in  the  engineering  literature.  In  addition  to  providing  a  solid  theoretical 
foundation  for  such  schemes,  we  have  systematically  tested  them  numerically  on  a 
number  of  examples  and  reported  some  of  our  findings.  One  result  of  these  investiga¬ 
tions  has  been  our  efforts  to  develop  alternate  schemes.  The  approximation  framework 
can  be  used  efficiently  to  develop  a  class  of  schemes  based  on  spline  or  “finite-element” 
ideas.  A  discussion  of  our  findings  in  this  regard  will  appear  in  a  manuscript  that  is 
currently  in  preparation. 

We  dose  with  several  further  remarks  that  we  have  added  in  the  final  version  of 
this  paper,  partly  in  response  to  referees’ queries  and  partly  as  a  result  of  our  subsequent 
efforts  and  findings  in  related  investigations.  First,  as  we  noted  in  Remark  4.1,  the 
generality  of  our  theoretical  framework  (q  dependent  spaces,  norms,  etc.)  is  not 
essential  to  treat  Example  4.1  or,  indeed,  any  of  the  specific  examples  discussed  above. 
However,  if  one  considers  a  parabolic  system  as  in  Example  4.2  for  which  the  function 
k  is  parameter  dependent,  the  q  dependence  of  the  appropriate  inner  product  is 
essential.  In  fact,  such  problems  arise  naturally  in  estimation  questions  for  porous 
madia  problems,  where  one  of  die  parameters  to  be  estimated  is  the  function  k  (the 
field  porosity)  itself.  A  treatment  using  the  theoretical  framework  developed  above 
in  connection  with  cubic  spline  approximations  is  outfitted  for  such  problems  in  [42]. 

With  regard  to  general  spline  approximation  schemes,  we  have,  since  this  paper 
was  first  written,  completed  certain  efforts  on  spfine-besed  techniques  (referred  to 
several  times  above)  in  the  context  of  the  theoretical  framework  given  above.  Second- 
order  parabolic  and  hyperbolic  systems  [43],  as  well  as  higher-order  equations  arising 
in  elasticity  [44],  have  been  treated  and  our  findings  have  been  most  positive  from 
both  computational  and  theoretical  viewpoints. 
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well,  was  also  helpful  in  various  discussions  on  the  practical  use  of  the  algorithms 
employed  in  connection  with  the  numerical  results  of  f  6, 
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